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Abstract 

We consider the two matrix model with an even quartic potential 
W{y) = 2/^/4 + ay'^/2 and an even polynomial potential V{x). The 
main result of the paper is the formulation of a vector equilibrium 
problem for the limiting mean density for the eigenvalues of one of 
the matrices Mi. The vector equilibrium problem is defined for three 
measures, with external fields on the first and third measures and an 
upper constraint on the second measure. The proof is based on a 
steepest descent analysis of a 4 x 4 matrix valued Riemann-Hilbert 
problem that characterizes the correlation kernel for the eigenvalues of 
Ml. Our results generalize earlier results for the case a = 0, where the 
external field on the third measure was not present. 
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1 Introduction and statement of results 
1.1 Hermitian two matrix model 

The Hermitian two-matrix model is a probability measure of the form 

(1.1) ^ exp (-nTr(F(Mi) + W{M2) - rMiMs)) dMidMa, 

defined on the space of pairs (Mi,M2) of n x n Hermitian matrices. The 



constant Z„ in (1.1) is a normalization constant, r G M\{0} is the coupling 



constant and dMidM2 is the flat Lebesgue measure on the space of pairs 



of Hermitian matrices. In (1.1), V and W are the potentials of the matrix 
model. In this paper, we assume ^ to be a general even polynomial and we 
take W to be the even quartic polynomial 

(1.2) w{y) = -^y^ + ^y\ a G M. 

Without loss of generality we may (and do) assume that 

(1.3) r>0. 

We are interested in describing the eigenvalues of Mi in the large n limit. 

In 05] the case a = was studied in detail. An important ingredient 
in the analysis of [H] was a vector equilibrium problem that describes the 
limiting mean eigenvalue distribution of Mi. In this paper we extend the 
vector equilibrium problem to the case a 7^ 0. 
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1.2 Background 



The two-matrix model (1.1) with polynomial potentials V and W was in- 
troduced in [591 EO] as a model for quantum gravity and string theory. The 
interest is in the double scaling limit for critical potentials. It is generally 
believed that the two-matrix model is able to describe all {p, q) conformal 
minimal models, whereas the one-matrix model is limited to {p, 2) minimal 
models \30\ BT| 148] . In [61] the two- matrix model was proposed for the 
study of the Ising model on a random surface, where the logarithm of the 
partition function (i.e., the normalizing constant Zn in ( jl.l] )) is expected to 
be the generating function in the enumeration of graphs on surfaces. For 
more information and background on the physical interest we refer to the 
the surveys [391 SO] , and more recent physical papers [9l HH [511 [52] 

The two matrix model have a very rich integrable structure that is con- 
nected to biorthogonal polynomials, isomonodromy deformations, Riemann- 
Hilbert problems and integrable equations, see e.g. [2] [TO ] [T2l [T3l ITU H6l [50 1 
[60l [66] . This is the basis of the mathematical treatment of the two matrix 
model, see also the survey [Tl]. 

The eigenvalues of the matrices Mi and M2 in the two-matrix model are 
a determinantal point process with correlation kernels that are expressed in 
terms of biorthogonal polynomials. These are two families of monic poly- 
nomials {pk,ni^)}'^=o and {qi,niy)}'^=o, where 

Pk,n has degree k and qi n has 

degree satisfying the condition 

(1.4) / / PM(^)9^,n(y)e-"(^(^)+^(^)-™^)dxd2/ = hljk,i. 



—00 J —00 



The polynomials are well-defined by (1.4) and have simple and real zeros 



]. Moreover, the zeros of pk,n and pk+i^n, and those of qi^n and qi+i,n are 
interlacing [43] . 

The kernels are expressed in terms of these biorthogonal polynomials 
and their transformed functions 

/oo 
'?/,n(y)e-"(^(^)-™^)^iy, 
-00 

/oo 
-00 

as follows: 

n— 1 ^ 

(1.5) K['1\xi,X2) = j^Pk,n{^l)Qk,n{^2), 

k=0 
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rt— 1 ^ 

;i.6) K[2\x,y) = -7o-Pk,n{x)qk,n{y), 



hi 

k=0 k,n 
rt-1 



[1.7) K^^{y,x) = -J^Pk,n{y)Qk,n{ 



-n{V (x)+W {y)—Txy) 



k=0 k,n 
rt-1 



;i.8) K^f{yi,y2) = ^ -T^Pk,n{yi)qk,n{V: 



hi 

k=0 k,n 

Then Eynard and Mehta [50l |72], see also [211 EZl E], showed that the 
joint probabihty density function for the eigenvalues xi, . . . , x„ of Mi and 
yi,. . . ,yn of M2 is given by 



Vixi, ...,Xn,yi,.. ■,yn) 



(n!)2 



det 



^ ^ K^^\xi, Xj) 



and the marginal densities take the form 



Ku\xi,yj) 



(^K^'l\yi,Xj)^_^^ (^K^22{yi,yj] 



=1/ 



:i.9) 



Vixi, . . . ,Xn,yi, ■ ■ . , yn)dxk+i ■ ■ ■ dxn %+i ■ ■ ■ dyn 



n—k+n—l times 

_ {n-k)\{n-l)\ 



det 



I^ll {Xi^ Xj ^ 
^21 -^jO 



kA 



?J=1 



In particular, by taking / = 0, so that we average over the eigenvalues 
yi, . . . ,yn of M2 , we find that the eigenvalues of Mi are a determinantal 



(rt) 

point process with kernel -f^ii , see (1.5). This kernel is constructed out of 
the biorthogonal family {pk,n}^=o and {Qi,n}^o and the associated deter- 
minantal point process is an example of a biorthogonal ensemble [23] . It is 
also an example of a multiple orthogonal polynomial ensemble in the sense 
of [62]. 

In order to describe the behavior of the eigenvalues in the large n limit, 
one needs to control the kernels (1.5)-(1.8) as n — t- 00. Due to special 



recurrence relations satisfied by the biorthogonal polynomials, there exist 



Christoffel-Darboux-type formulas that express the n-term sums (1.5) and 



(1.8) into a finite number (independent of n) of biorthogonal polynomials 
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and transformed functions, see [13]. This paper also gives differential equa- 
tions and a remarkable duality between spectral curves, see also [T^ 

A Riemann-Hilbert problem for biorthogonal polynomials was first for- 
mulated in [46\ . The Riemann-Hilbert problem in [46j is of size 2x2 but it 
is non-local and one has not been able to apply an asymptotic analysis to it. 
Local Riemann-Hilbert problems were formulated in |14^ l60l l66] , but these 
Riemann-Hilbert problems are of larger size, depending on the degrees of the 
potentials V and W. The formulation of a local Riemann-Hilbert problem 
for biorthogonal polynomials, however, opens up the way for the application 
of the Deift-Zhou [33] steepest descent method, which was applied very suc- 
cessfully to the Riemann-Hilbert problem for orthogonal polynomials, see 
[THl [5H 133 [Mj and many later papers. 

In [35] the Deift-Zhou steepest descent method was indeed applied to 
the Riemann-Hilbert problem from [66j for the case where W is given by 
(1.2) with a = 0. It gave a precise asymptotic analysis of the kernel 
as n — 7- oo, leading in particular to the local universality results that are 
well-known in one-matrix models j33j . The analysis in [35] was restricted to 
the genus zero case. The extension to higher genus was done in [75j . 



1.3 Vector equilibrium problem 

As already stated, it is the purpose of the present paper to extend the results 
of [m [75] to the case of general a. 

An important role in the analysis in |45j is played by a vector equilibrium 
problem that characterizes the limiting mean density for the eigenvalues 
of Ml (and also gives the limiting zero distribution of the biorthogonal 
polynomials Pn,n)- One of the main contributions of the present paper is the 
formulation of the appropriate generalization to general a G M. We refer to 
the standard reference |78j for notions of logarithmic potential theory and 
equilibrium problems with external fields. 



1.3.1 Case a = 

Let us first recall the vector equilibrium problem from |45j . which involves 
the minimization of an energy functional over three measures. For a measure 
/i on C we define the logarithmic energy 

HfJ-) = / / log 1 — ^—rdn{x)dfi{y) 
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and for two measures fj, and u we define the mutual energy 

HfJ':^) = I I log] — - — \dii{x)dv{y). 
J J \x-y\ 

The energy functional in [33] then takes the form 



3 2 

:i.io) / {y{x)-\\Tx\''') 

i=i .7=1 ^ ^ 



du\[x) 



and the vector equilibrium problem is to minimize ( 1.10 ) among all measures 
1^1, V2 and f3 such that 



(a) the measures have finite logarithmic energy; 

(b) v\ is a measure on M with z^i(M) = 1; 

(c) V2 is a measure on iM with 1^2(^1^) = 2/3; 



(d) 1/3 is a measure on M with z^3(M) = 1/3; 



(e) V2 < o"2 where (T2 is the unbounded measure with density 

(1.11) ^ = ^^4/3|^|l/3^ ^^.^ 

\az\ 2tt 



on the imaginary axis. 



A main feature of the vector equilibrium problem is that it involves an 
external field acting on the first measure as well as an upper constraint 
( |1.11[ ) acting on the second measure. Note that an upper constraint arises 
typically in the asymptotic analysis of discrete orthogonal polynomials, see 
e.g. [3 |22l II2I [68l [77j . The interaction between the measures in (1.10) is of 
the Nikishin type where consecutive measures attract each other, but there 
is no direct interaction between measures Vi and Vj \i \i — j\ > 2. The 
notion of a Nikishin system originated in works on Hermite-Pade rational 
approximation, see [H [56l ESI ES] . Vector equilibrium problems also played 
a role in the recent papers [H [151 [13 ESj that are related to random matrix 
theory and [Gj [131 HU [691 El] that are related to recurrence relations. 
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1.3.2 General a 

For general a G M, the relevant energy functional takes the form 

3 2 

(1.12) E{u^,iJ2,i^3) = Y,I{i^j) -Y,I{v,,v,^i) 

+ j Vi{x)dvi{x) + j Vi{x)dv3{x), 

where Vi and V3 are certain external fields acting on vi and z/3, respec- 
tively. The vector equilibrium problem is to minimize £{1^1, 1^2, 1^3) among 
all measures 1^1, 1^2, fs, such that 

(a) the measures have finite logarithmic energy; 

(b) vi is a measure on M with 1^1 (M) = 1; 

(c) 1^2 is a measure on iR with 1^2(^1^) = 2/3; 

(d) z^3 is a measure on M with 1^3(1^) = 1/3; 

(e) 1^2 ^ o"2 where (T2 is a certain measure on the imaginary axis. 



Comparing with (1.10) we see that there is an external field V3 acting 
on the third measure as well. The vector equilibrium problem depends on 
the input data Vi, V3, and (T2 that will be described next. Recall that V is 



an even polynomial and that W is the quartic polynomial given by (1.2) 



External field Vi: The external field Vi that acts on z^i is defined by 

(1.13) Vi(x) = V(x) + mm (W(s) - Txs) . 

The minimum is attained at a value s = si{x) G M for which W'{s) = tx, 
that is 

(1.14) s^ + as = Tx. 



For a > 0, this value of s is uniquely defined by (1.14). For a < there 



can be more than one real solution of (1.14). The relevant value is the one 
that has the same sign as x (since r > 0, see (1.3)). It is uniquely defined, 
except for X = 0. 
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External field V3: The external field V3 that acts on 1/3 is not present if 
Q > 0. Thus 



(1.15) 



= 



if a > 0. 



For a < 0, the external field V3(x) is non-zero only for x G {—x*{a),x*{a)) 
where 



:i.i6) 



x*{a) 



-a 



3/2 



a < 0, 
a > 0. 



For those x, the equation (1.14) has three real solutions si = si{x), S2 = 
S2{x), S3 = S3(x) which we take to be ordered such that 

W{si) — TXSl < W{S2) — TXS2 < W{s3) — TXS^. 

Thus the global minimum of s G M 1— )■ W{s) = txs is attained at si, and this 



global minimum played a role in the definition (1.13) of Vi. The function 



has another local minimum at S2 and a local maximum at S3, and these are 

used in the definition of V3. We define V3 : M — )■ M by 

(1.17) 

'(l^(s3(x)) - rxs3(3;)) - 

{W{s2{x)) — TXS2{x)) , for X E (— a;*(a), x*(a)), 
elsewhere. 



Thus V^{x) is the difference between the local maximum and the other local 
minimum of s G M 1— )■ W{s) = txs, which indeed exist if and only if x G 



{—x*{a),x*{a)), where x*{a) is given by (1.16). In particular Vs{x) > for 
X G (— x*(q), x*{a)). 



The constraint CJ2: To describe the measure (T2 that acts as a constraint 
on 1^2, we consider the equation 



:i.i8) 



+ as = Tz, with z G iH 



There is always a solution s on the imaginary axis. The other two solutions 
are either on the imaginary axis as well, or they are off the imaginary axis, 
and lie symmetrically with respect to the imaginary axis. We define 



(1.19) 



daoiz) T 



\dz\ 



vr 



Re s{z) 
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where s{z) is the solution of (1.18) with largest real part. We then have for 
the support S{(J2) of c72, 

(1.20) S{a2)=iR\{-iy*{a),iy*{a)), 
where 

r 2 /a\3/2 

(1.21) y*{a) = 



a > 0, 



T V3. 
0, a < 0. 



This completes the description of the vector equilibrium problem for 
general a. It is easy to check that for a = it reduces to the vector 
equilibrium described before. 

1.4 Solution of vector equilibrium problem 

Our first main theorem deals with the solution of the vector equilibrium 
problem. We use S{fi) to denote the support of a measure The logarith- 
mic potential of /x is the function 

(1.22) Wix) = / log , ^ , dfi{s), xeC, 



\x — s\ 

which is a harmonic function on C \ 5(^) and superharmonic on C 

Theorem 1.1. The above vector equilibrium problem has a unique mini- 
mizer (^i,/X2,/i3) that satisfies the following. 

(a) There is a constant S M such that 
(1.23) 

//O S{fii) or 5(c72 - 1^2) then 



2U>''{x) = U^'^ix) -Vi{x) + £i, xe S{fii), 
2U^^{x) > U''^{x)-Vi{x) + ii, x£R\S{ni). 



N 

(1.24) S{i,i) = [J[a„bj], 

for some N £ N and ai < bi < 02 < • • ■ < on < bjy, and on each of 
the intervals [aj,bj] in S{fii) there is a density 



(1-25) ^ = pi{x) = ^hj{x)^J{bj -x){x-aj), x G [aj, 
and hj is non-negative and real analytic on [aj,bj]. 
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and there is a constant C2 > such that 
(1.27) Sifi2) = S{a2), and 5(a2 
Moreover, a2 — has a density 
d{(J2 - ^2) 



(1.28) 



\dz\ 



P2{z), 



X G S{a2 - /i2), 

X eiM. \ S{a2 - ^2), 



fi2) = iM \ {-ic2, ic2). 



z G 



that is positive and real analytic on iM \ [—ic2,ic2]- If C2 > 0, then 
P2 vanishes as a square root at z = ibic2. If a>0, then C2 > y*{a), 
where y*{a) is given by (1.21). 

(c) We have 



(1.29) 



2U^''{x) = U^'^ix) -Vsix), 
2U^^{x) > U^'^x) -V3{x), 



X e S{fi3), 
X eR\S{fis), 



and there is a constant C3 > such that 

(1.30) 5(/i3) =M\(-C3,C3). 

Moreover, 113 has a density 

(1.31) ^ = P3ix), 



X G 



that is positive and real analytic on ]R\ [— C3, C3]. If a >0, then C3 = 0. 
If a < 0, then C3 < x*{a) where x*{a) is given by (1.16). // C3 > 0, 
then p3 vanishes as a square root at x = ±03. 

(d) All three measures are symmetric with respect to 0, so that for j = 
1, 2, 3 we have f^j{A) = fij{—A) for every Borel set A. 

In part (a) of the theorem it is stated that S{fii) is a finite union of 
intervals under the condition that S{fj,i) and S{a2 — P2) are disjoint. If this 
condition is not satisfied then we are in one of the singular cases that will 
be discussed in Section 1.5 below. However, the condition is not necessary 
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as will be explained in Remark 4.9 below. We chose to include the condition 



in Theorem 1.1 since the focus of the present paper is on the regular cases. 



The conditions (1.23), (1.26), and (1.29) are the Euler-Lagrange varia- 



tional conditions associated with the vector equilibrium problem. We note 



the strict inequalities in (1.26) and (1.29). These are consequences of special 



properties of the constraint a2 and the external field V3 that are listed in 
parts (b) and (c) of the following lemma. 

Lemma 1.2. The following hold. 

(a) Let U2 he a measure on iM such that V2 < (J2- If ^ ^ S{(T2 — 1^2) then 
X I— )• Vi{x) — ?7^2(x) is real analytic on R. 



(b) The density f^Jiy) = ^Resiiy) (see (L19\) ) is an increasing function 



for y > 0. 



(c) Let a < 0. Let U2 he a measure on iM of finite logarithmic energy such 
that 1^2 < (^2- Then x ^ V'i{y/x) — U^'^{y/x) is a decreasing and convex 
function on (0, (x*(q))^). 



IS given 



Lemma 1.2 is proved in Section [2] and the proof of Theorem 1.1 
in Section [3l 

A major role in what follows will be played by functions defined on a 
compact four-sheeted Riemann surface that we will introduce in Section |4j 
The sheets are connected along the supports S{^i), S{a2 — ^2) and S{iJ,3) 
of the minimizing measures for the vector equilibrium problem. The main 
result of Section |4] is Proposition 4.8 which says that the function defined 
by 

V'{z) - [ ^ 
J z - s 

on the first sheet has an extension to a globally meromorphic function on 
the full Riemann surface. This very special property is due to the special 
forms of the external fields Vi and V3 and the constraint a2 , which interact 
in a very precise way. 



1.5 Classification into cases 

According to Theorem |l.l| the structure of the supports is the same for a > 
as it was for a = in [45J, that is, S{fj,3) = M and S'(cr2— /X2) = ZC2, ^02) 

for some C2 > 0. The supports determine the underlying Riemann surface, 
and so the case a > is very similar to the case a = 0. There are no 
phase transitions in case a > 0, except for the possible closing or opening 
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of gaps in the support of //i. These type of transitions already occur in the 
one- matrix model. 

For a < 0, however, certain new phenomena occur which come from the 
fact that the external field V3 on (defined in (1.17)) has its maximum at 
and therefore tends to move fi^ away from 0. As a result there are cases 
where S^fi^) is no longer the full real axis, but a strict subset (1.30) with 
C3 > 0. 

In addition, it is also possible that C2 = in (1.27) such that 5(ct2 — ^2) 
is the full imaginary axis and the constraint (T2 is not active. These new 
phenomena already occur for the simplest case 

Vix) = 

for which explicit calculations were done in [53] based on the coefficients in 
the recurrence relations satisfied by the biorthogonal polynomials. These 
calculations lead to the phase diagram shown in Figure |1.1| which is taken 
from [33]. There are four phases corresponding to the following four cases 
that are determined by the fact whether is in the support of the measures 
//I, fT2 - ^12, Ms or not: 

Case I: G ^(/xi), S{a2 - 1^2), and G 3(^13), 

Case II: S{a2 - ^2), and G ^(/ia), 

Case III: S{ni), G S{a2 - ^2), and 5(^3), 

Case IV: G S{fii), ^ S{a2 - /i2), and 8(113). 

The four cases correspond to regular behavior of the supports at 0. There 
is another regular situation (which does not occur for V{x) = jx^), namely 

Case V: S{fii), S{a2 - /X2), and ^(/is). 

The five cases determine the cut structure of the Riemann surface and we 
will use this classification throughout the paper. 

Singular behavior occurs when two consecutive supports intersect at 0. 

Singular supports I: G S{fii) Pi S{a2 — IJ-2), S^Js), 

Singular supports II: S{fii), G S{a2 - IJ.2) n S{fi3). 

There is a multisingular case, when all three supports meet at 0: 



Singular supports III: G ^(^i) n S{a2 - /X2) n ^(/is). 
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Figure 1.1: Phase diagram in the a-r plane for the case V{x) = ^x^: the 
curves r = ^/a+~2 and r = \J —\ja separate the phase diagram into four 
regions. The four regions correspond to the cases: Case I: = 1, C2 > 
and C3 = 0, Case II: iV = 2, C2 > and C3 = 0, Case III: = 2, C2 = and 
C3 > 0, and Case IV: = 1, C2 > and C3 > 0. 



Besides a singular cut structure for the Riemann surface, we can also 
have a singular behavior of the first measure These singular cases also 
appear in the usual equilibrium problem for the one-matrix model, see [33] , 
and they are as follows. 

Singular interior point for The density of \x\ vanishes at an interior 
point of S{fii). 

Singular endpoint for /ii: The density of /zi vanishes to higher order 
than square root at an endpoint of S{p,i). 

Singular exterior point for Equality holds in the variational inequal- 



ity in ( 1.23 ) at a point 2; G M \ S{fii). 



The measures cj2 — /i2 and /i3 cannot have singular endpoints, singular 
exterior points, or singular interior points, except at 0. Singular interior 
points of these measures at are as follows. 



Singular interior point for a2 — fJ'2' The density of ct2 — /i2 vanishes at 
G S{a2 — /U2). 
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Singular interior point for The density of ^3 vanishes at G ^(/Us). 

While there is great interest in the singular cases we restrict the analysis 
in this paper to the regular cases, for which we give the following precise 
definition. 

Definition 1.3. The triplet {V,W,t) is regular if the supports of the mini- 
mizers from the vector equilibrium problem satisfy 

(1.32) 5(^1) n 5(^72 -;U2) =0 and ^(/is) n 5(cj2 - /U2) = 

and if in addition, the measure /Ui has no singular interior points, singular 
endpoints, or singular exterior points, and the measures (T2 — ^2 'I'^-d fJ-s do 
not have a singular interior point at 0. 



The condition (1.32) may be reformulated as 



C2 = =^ O0S(^i)U5(^3). 

1.6 Limiting mean eigenvalue distribution 

The measure /Ui is the limiting mean eigenvalue distribution of the matrix 
Ml in the two- matrix model as n — )• 00. In this paper we prove this only for 
regular cases. To prove it for singular cases, one would have to analyze the 
nature of the singular behavior which is beyond the scope of what we want 
to do in this paper. 

Theorem 1.4. Suppose {V,W,t) is regular. Let /xi be the first component 
of the minimizer [m, 112, fj-s) of the vector equilibrium problem. Then fii 
is the limiting mean distribution of the eigenvalues of Mi as n — )• 00 with 
n = 0(mod3). 

We recall that the eigenvalues of Mi after averaging over M2 are a de- 



terminantal point process on M with a kernel i^j^ as given in (1.5). The 
statement of Theorem 11.41 comes down to the statement that 



(1.33) lim ^k\'^\x,x) = pi(x), x€R 

n— >oo n 

where pi is the density of the measure pi. 

The restriction to n = 0(mod3) is for convenience only, since it simpli- 
fies the expressions in the steepest descent analysis of the Riemann-Hilbert 
problem that we are going to do. 



16 



M. Duits, A.B.J. Kuijlaars, and M.Y. Mo 



The existence of the Hmiting mean eigenvalue distribution was proved 
by Guionnet [S^ in much more general context. She characterized the min- 
imizer by a completely different variational problem, and also connects it 
with a large deviation principle. It would be very interesting to see the 
connection with our vector equilibrium problem. A related question would 
be to ask if it is possible to establish a large deviation principle with the 



energy functional (1.12) as a good rate function. 



We are going to prove ( 1.33 ) by applying the Deift-Zhou steepest descent 



analysis to the Riemann-Hilbert problem (1.36) below. Without too much 



extra effort we can also obtain the usual universal local scaling limits that 
are typical for unitary random matrix ensembles. Namely, if pi{x*) > 
then the scaling limit is the sine kernel 



1 



lim 

n^oo npi{x*) 



,x + 



y 



smvr X 



y) 



npi{x*)^ ' npi{x*) ) ■k{x — y) 

while if X* G {oi, 6i, . . . , oat, 6Ar} is an end point of S{pi) then the scaling 
is the Airy kernel, i.e., for some c > 0, we have 

1 .Jn) / . . X ^ . y \ Ai(x) Ai'(y) - Ai'(x)Ai(y) 



lim 



:K['1^ ( x" ± 



=P I 

, X ± 



n->oo (^cn)^/'^ \ (cn)-^/'^' (cn)2/^y x — y 

with + if X* = bj and — if x* = aj for some j = 1, . . . ,N. Recall that we 
are in the regular case so that the density of pi vanishes as a square root at 
x* . The proofs of these local scaling limits will be omitted here, as they are 
very similar to the proofs in pBj. 



1.7 About the proof of Theorem 1.4 



The first step in the proof of Theorem 1.4 is the setup of the Riemann- 
Hilbert (RH) problem for biorthogonal polynomial pn,n and its connection 
with the correlation kernel k[^\ We use the RH problem of |66j which we 
briefly recall. 

The RH problem of |66j is based on the observation that the polynomial 
Pnn that is characterized by the biorthogonality conditions (1.4) can alter- 



natively be characterized by the conditions (we assume W is quartic and n 
is a multiple of three) 



(1.34) 



Pn,n{x)x^Wj^n{x)dx = 0, /c = 0, . . . , n/3 - 1, 3 = 0, 1, 2, 



which involves three varying (i.e., n-dependent) weight functions 



(1.35) 



w 



X 



yj Q-n{W{y)-Txy) 



dy, 



0,1,2. 
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The conditions (1.34) are known as multiple orthogonality conditions of type 
II, see e.g. [H ESI 113 El- 

A RH problem for multiple orthogonal polynomials was given by Van 
Assche, Geronimo and Kuijlaars in as an extension of the well-known 
RH problem for orthogonal polynomials of Fokas, Its, and Kitaev [55]. For 
the multiple orthogonality (1.34) the RH problem is of size 4x4 and it asks 
for y : C \ M ^ C^^^ satisfying 



(1.36) < 



y+(x) = y„(x) 



Y{z) = {I + Oil/z)) 



Y is analytic in C \ M, 

/I Wo,n{x) Wi^n{x) W2,nix)\ 

1 

1 

Vo 1 / 

/z" \ 




Vo z-"/V 



X G 



z — )• oo. 



The RH problem has a unique solution. The first row of Y is given in 
terms of the biorthogonal polynomial pn,n as follows 



Yl,l{z) = Pn,niz), Yij+2{z) 



1 

2TTi 



Pn,n{x)wj^n{x) 



X — z 



dx, j = 0,1,2, 



and the other rows are built out of certain polynomials of degree n — 1 in a 
similar way, see [66l ES] for details. 

Multiple orthogonal polynomials have a Christoffel-Darboux formula 



which implies that the correlation kernel (1.5) can be rewritten in the inte- 
grable form 

fi{x)gi{y) + h{x)g2{y) + h{x)g^{y) + fi{x)gA{v) 
x-y 

for certain functions fj,gj, for j = 1, . . . , 4, and in fact it has the following 
representation 



(1.37) K['l>{x,y) 



1 

2TTi{x — y) 



(o wo,n(y) wi,n{y) w2,n{y))^+^iy)Y+(.x) 



n\ 
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for x,y gM, in terms of the solution Y of the RH problem (1.36), see |29j . 

The proof of Theorem 1.4 is an involved and lengthy steepest descent 
analysis of the RH problem 1.36 in which the vector equilibrium problem is 
used in an essential way. This is similar to |45j which deals with the case 
a = 0. Certain complications arise because the formulas for the external 
field and the constraint in the vector equilibrium problem are less explicit 
as in the case a = 0. This not only complicates the analysis of the vector 
equilibrium problem in Sections [2] and [3} but it will continue to play a role 
via the functions 6j defined in Section |2.2| and A 
throughout the paper. 

We also note that the analysis in 



defined in Section 4.4 



was restricted to the one-cut case, 
which leads to an underlying Riemann surface of genus 0. This restriction 
was removed in [75] . The problem in the higher genus case is in the construc- 
tion of the global parametrix. In Section [8] we give a self-contained account 
that is based on the ideas developed in j75] and [67], which we think is of 
independent interest. 

We also wish to stress that in Case IV the Riemann surface always has 
genus > 1, even if S{fii) consists of one interval, see (4.16) below. This 
phenomenon did not appear for q = 0. 



1.8 Singular cases 

Although we do not treat the singular cases in this paper we wish to make 
some comments about the possible critical behaviors that we see in the two- 
matrix model with the quartic potential W{y) = + ^y^. 

As already discussed in Section |1.5| the singular behavior is associated 
with either a singular behavior in the measures //i, ct2 — IJ^2, or //s, or a 
singular behavior in the supports. The singular behavior in the measure 
Hi also appears in the one-matrix model that is described by orthogonal 
polynomials. It is known that the critical behavior at a singular interior 
point where the density vanishes quadratically is described by the Hastings- 
McLeod solution of the Painleve II equation, see |19 ^ 1271 179j . This Painleve 
II transition is the canonical mechanism by which a gap opens up in the 
support in the one-matrix model. 

The critical behavior at a singular endpoint where the density vanishes 
with exponent 5/2 is described by a special solution of the Painleve I2 equa- 
tion (the second member of the Painleve I hierarchy), see [28]. The criti- 
cal behavior at a singular exterior point is described by Hermite functions 

126] [73] and this describes an opening of a new band of eigenvalues (birth 
of a cut). 
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We see these critical behaviors also in the two-matrix model with an 
even quartic W. In particular, the opening of a gap at in the support of 
fii is a Painleve II transition. In our classification of regular cases, this is a 
transition from Case I to Case II, or a transition from Case IV to Case V. 



In the phase diagram of Figure 1.1 for V{y) = ^y'^, this transition is on the 



part of the parabola r = \/a + 2, with a > —1. 

A Painleve II transition also appears when either CJ2 — //2 or ^3 has a 
density that vanishes quadratically at 0. Then is a singular interior point 
and again a gap can open but now in the support of the measures "that are 
on the other sheets" and have no direct probabilistic meaning. If the density 
of £72 — //2 vanishes at then the transition is from Case III to Case V. If 
the density of /Us vanishes at then the transition is from Case I to Case IV 
or from Case II to Case V. In the phase diagram of Figure [LT] the transition 
from Case I to Case IV takes place on the part of the parabola r = \/a + 2, 
with -2 < a < -1. 

The cases of singular supports represent critical phenomena that do not 
appear in the one-matrix model. What we called Singular Supports I in 
Section 1.5 corresponds to a transition from Case III to Case IV. This is 
a transition when the gap around in the support of S{ni) closes and 
simultaneously the gap in the support of S{a2 — ^12) opens up (or vice versa). 
On the level of the Riemann surface it means that the two branch points on 
the real line that are the endpoints of the gap in S{fii) come together at 0, 
and then split again to become a pair of complex conjugate branch points. 
These branch points are then on the imaginary axis and are the endpoints 
±ic2 of S{a2 — ^2)- A transition of this type does not change the genus of 
the Riemann surface. 

This type of transition was observed first in the context of random 
matrices with external source and non-intersecting Brownian motions, see 
O EU EH [8l],i where it was described in terms of Pearcey integrals. The 
Pearcey transition is a second mechanism by which a gap in the support 
may open up (or close). As it involves three sheets of the Riemann sur- 
face it cannot take place in the one-matrix model which is connected to a 
two-sheeted Riemann surface. 

The case of Singular Supports II gives a transition from Case II to Case 
III. This is a situation where the gap in the support of (T2 — H2 closes and 
simultaneously the gap in 8(113) opens. This also typically corresponds to 
a Pearcey transition, but it does not involve the first sheet of the Riemann 
surface, which means that this transition is not visible in the eigenvalue 
distribution of the random matrix. In the phase diagram of Figure |1.1| the 
Pearcey transitions are on the curve r = y^—l/a, a ^ —1. 
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The case of Singular Supports III represents a new critical phenomenon. 
Here the supports of all three measures cr2 — ^2 and fi^ are closed at 0. In 



Figure 1.1 this is the case at the multi-critical point a = —1 and r = 1 where 
the Painleve transitions and Pearcey transitions come together. One may 
approaches the multi-critical point from the Case III region, where there 
is a gap around in the supports of both fii and /is, while the support of 
^2 — fJ'2 is the full imaginary axis. At the multi-critical point the supports 
of fii and H3 close simultaneously, while also the support of a"2 — H2 opens 
up, which results in a transition from Case III to Case I. 

We conjecture that the case of Singular Supports III is of similar nature 
as was studied very recently [T, "38] for a critcal case of non-intersecting 
Brownian motions (or random walks) with two starting and two ending 
points. By fine-tuning the starting and ending points one may create a 
situation where two groups of non-intersecting Brownian motions fill out 
two ellipses which are tangent to each other at one point. Our conjecture 
is that the local eigenvalue behavior around in the multi-critical case is 
the same as that for the non-intersecting Brownian motions at the point 
of tangency. The conjecture is supported by preliminary calculations that 
suggest that the local parametrix of [38] can also be used if one tries to 
extend the RH analysis of the present paper to the multi-critical situation. 



2 Preliminaries and the proof of Lemma 1.2 



Before coming to the proof of Theorem |1.1| we study the equation (1.14) in 



more detail. This equation will also play a role in the proof of Theorem 1.4 
where in the first step of the steepest descent analysis, we will use functions 
defined by integrals 

(2.1) / e-"('^('^)-^^^)ds, W{s) = + ^s^ 

Jr 4 2 

where T is an unbounded contour in the complex z-plane. 
2.1 Saddle point equation and functions Sj 



The large n asymptotics of the integrals (2.1 ) is determined by the solutions 



of the saddle point equation W'{s) — tz = 0, that is 



(2.2) 



+ as = TZ. 
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In (1.14) we considered this equation for z = x £ M. We defined a solution 
si{x) for every x G M, and for a < and |x| < x*{a) we also defined S2{x) 
and S3(x). 



We define solution si{z), S2{z) and 53(2;) of (2.2) for complex z as follows. 
We distinguish between the two cases a > and a < 0. 



Case a > 0. In case a > the saddle point equation (2.2) has branch 



points ziziy*{a) G iM where y*{a) is given by (1.21). The Riemann surface 



S for the equation (2.2) then has three sheets that we choose as follows 



(2.3) 



' cSi = C \ {{-ioo, -iy*{a)] U [iy* {a) , ioo)) , 

52 = C \ (M U (-ioo, -iy*{a)] U [iy*(a), ioo)) , 

53 = C \ M. 



We already defined si{x) for x G M as the unique real saddle point. This 
function has an analytic continuation to Si that we also denote by si. Then 
S2 and S3 are defined by analytic continuation onto ^2 and ^3, respectively. 



Case a < 0. In case a < the saddle point equation (2.2) has two real 



branch points ±x*(a) with x*{a) given by (1.16). The three sheets of the 



Riemann surface S for the equation (2.2) are now chosen as follows 



(2.4) 



52 = C \ ((-00, -x*{a)] U [x*{a),oo) U iR) , 

53 = C \ ((-00, -x*{a)] U [x*{a),oo)) . 



In case a < 0, we have that si{x) is defined for x G M \ {0}. It is the 
real saddle point for which W{s) — txs is minimal. The function si has 
an analytic continuation to Si that we also denote by si. Then S2 and S3 
are defined by analytic continuation onto ^2 and ^3, respectively. It is a 
straightforward check that for x G {—x* (a), x* (a)) this definition of S2(x) 
and S3(x) coincides with the one earlier given. 

Lemma 2.1. The functions Sj are have the symmetries 



(2.5) Sj{-z) = -Sj{z), Sj{z) = Sj{z), j = 1,2,3. 

In addition we have that 



(2.6) 



Resi(z) > if Rez > 0. 
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Proof. The symmetries (2.5) are clear. 

For z = X £ M. with x > we have that si{x) > 0. Therefore, by 
continuity, Jiesi{z) > for z in a neighborhood of the positive real axis. If 
Resi(z) = for some z, so that si{z) is purely imaginary, then 

TZ = si{z)'^ + asi{z) 

is purely imaginary as well. The inequality Re ,81(2:) > therefore extends 
into the full right half-plane as claimed in (2.6). □ 



From the lemma it follows that in both cases the constraint cr2, see (1.19), 
is given by 

dcr2{z) T 

— = — Resi -(z) 

az TTl 

(2.7) = ^(si,_(z)-si,+ (z)), zGiR. 

The imaginary axis is oriented upwards, so that si-{z) (si,-^-{z)) for z G 
denotes the limiting value of si as we approach z £ iR from the right (left) 
half-plane. 



2.2 Values at the saddles and functions 6', 



We define 

(2.8) 9j{z) = -W{sj{z))+Tz.Sj{z), j = 1,2,3 

as the value of —{W{.s) — tzs) at the saddle s = Sj{z). Note that 

(2.9) 9'j{z) = {-W'isjiz)) + Tz) s'jiz) + TSj{z) = rsjiz) 

so that, up to a factor r, 9j is a primitive function of sj. 

Then 9j is defined and analytic on Sj, see (2.3) and (2.4), and 



(2.10) 



92,zp on {—ioo, —iy*{a)] U [iy*{a), ioo), 
93^=p on (— 00, — x*(a)] U [x*(a), 00). 



Recall from (1.16) and (1.21) that we have put x*{a) = if a > and 
y*{a) = if a < 0, so that we can treat the two cases simultaneously in 

(|27To|. 
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The jumps for the 6j functions from (2.10), are taken together in terms 
of the jumps of the diagonal matrix 



h{z) 
(2.11) G(z) = ( e2{z) 

Osiz), 

as follows. 

Corollary 2.2. For x gM. we have 



z G C \ (M U iM) 



(2.12) 



f G+(x) = e_(x), 

1 0\ /I 0> 

G+(x) = I 1 e_(x) 1 



^0 1 Oj 

For z = iy G iM we have 



.0 1 0> 



(2.13) 



e+(z) = G_(z), 

/O 1 0\ /O 1 0> 

e+(z) =1 @-{z) 10 

\0 1/ \0 1, 



|x| < x*(a), 
\x\ > x*{a). 

\y\ < y*{a), 
\y\ > y*{a). 



Also note that by (1.13), the definition of si and (2.8), we have 
(2.14) 



Vi{x) = V{x)-9i{x), 



X £ 



and by ( 1.15|)-(|1.17|), the definition of S2 and S3, and (|2.8| 
(2.15) 



2{x)-93{x), foT X e {-x*{a),x*{a)), 
0, elsewhere. 



2.3 Large z asymptotics 

In what follows we will need the behavior of Sj{z) and 9j{z) as 2 — )• 00. 

Throughout the paper we define fractional exponents with a branch cut 
along the negative real axis. We use /, //, /// and IV to denote the four 
quadrants of the complex z-plane. We also put 



u = e 



2ni/3 



We state the following lemma without proof. It follows easily from the 
saddle point equation (2.2). 
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oo 



si{z 



Lemma 2.3. We have as z 

'(tz)1/3 _ |(rz)-i/3 + |J(tz)-5/3 + 0(z-7/3), 

L^2(rz)l/3 _ |a;(r2)-l/3 + |!^2(^^)-5/3 + 0(^-7/3)^ 
L<j(tz)1/3 _ |^2(^^)-l/3 + |i^(^^)-5/3 + 0(^-7/3)^ 

(tz)1/3 - f (rz)-i/3 + |J(tz)-5/3 + 0(z-7/3), 

w2(rz)l/3 _ |^(^2)-l/3 ^ |!^2(^^)-5/3 ^ 0(^-^/3), 
a;2(^^)l/3 _ |^(^^)-l/3 + |i^2(^^)-5/3 + 0(^-7/3)^ 
a;(Tz)V3 _ |^2(^^)-l/3 + |i^(^^)-5/3 + 0(^-7/3)^ 



m I U IV, 
in II, 
in III, 

in I 

in II U ///, 
in IV, 

in lull 
in III U IV. 



We have a similar result for the asymptotics of 6j . Note that the following 
asymptotic behaviors are consistent with the property that 9j = tsj, see 



Lemma 2.4. We have 



as z 



oo 





fM2/3 + ^- 


54 


V^)-2/3 + 0(^-4/3)^ 




in 


luiv 


< |^^(rz)4/3- 


- f a;2(rz)2/3 + 


a2 
6 


- |it^(rz)-2/3 + 


4/3)^ 


in 


II, 


J^2(^,)4/3 


- f w(rz)2/3 + 


a2 
6 


- fJw2(Tz)-2/3 + 


-4/3)^ 


in 


III, 


'|^^(rz)4/3- 


- f a;2(rz)2/3 + 


a2 
6 


- gw(rz)-2/3 + 


4/3)^ 


in 


I 


< |(rz)V3_ 


fM2/3 + ^f - 


.-I 
54 


;rz)~2/3 + 0(^-4/3)^ 




in 


II u ///, 


J^2(^,)4/3 


- f 6^(rz)2/3 + 


a2 
6 


- fJw2(Tz)-2/3 + 0(z 


-4/3)^ 


in 


IV, 


f f^^2(^^)4/3 


- f a;(Tz)2/3 + 


a2 
6 


fjL^2(^^)-2/3 ^ 


-4/3)^ 


in 


lull 


\|^(rz)V3_ 


- f a;2(Tz)2/3 + 


a2 
6 


g^(rz)-2/3 + 0(z- 


4/3)^ 


in 


III U IV. 



2.4 Two special integrals 



As a final preparation for the proof of Lemma 1.2 we need the evaluation of 
the following two definite integrals. 



Lemma 2.5. We have for x > 
da2{z) 



(2.16) 



(x — z)^ 



-rs'i(x), 



and 



d(72{z) _ TSl{^/x) 



X — z^ 



Proof. Because of the formula (2.7) for (T2 we have 
da2{z) 



{x - z)"' 



T 
2TTi 



si,^{z) - Sl,+ (^) 
(x — z)^ 



dz 
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Since si is analytic in C\iM and si{z) = 0{z^^^) as z — )• c«, see Lemma 2.3 
we can evaluate the integral using contour integration and residue calculus. 
It follows that 



r 
27ri 



{x 



dz 



-TSi [X 



and 



271-i 



dz = 0, 



and the first integral in (2.16) is proved. 



The second integral follows by a similar calculation, where we also use 
the fact that si is an odd function. □ 

2.5 Proof of Lemma 11.21 



Now we come to the proof of Lemma 1.2 



2.5.1 Proof of part (a) 



Proof. Integrating the first formula in (2.16) two times with respect to x, 



and using the fact that 9'i = tsi, we find that for some constants A and B, 
/ {log\z — x\ —log\z\) da2{z) = 6i{x) + Ax + B, x > 0. 

JiM. 

Thus 



Vi{x) - U'^^ix) = V{x) - ei{x) + J log \z - x\ dv2{x) 

(2.17) =V{x)- I {\og\z-x\-\og\z\)d{a2-V2){z)-Ax-B' 

JiR 

with a different constant B' = B — J log \ z\ du2{z)- 

Since S{a2 — 1^2) there exists a c > such that S{a2 — 1^2) C «M \ 



{—ic,ic) and so the integral in the right-hand side of (2.17) defines a real 



analytic function of x. Then (2.17) proves that Vi — U^^ is real analytic on 
M, since y is a polynomial. □ 



2.5.2 Proof of part (b) 

Proof. We have by ([2jJ) 

da2 

W\ 



(iy) = -Resi_(zy) 



vr 
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Since si is a solution of + as = tz, we have (3sf + a)s'i = r, so that 



(2.18) 



Re si _ {iy) = -r Im ^ r \ , — • 



For y > we have Resi^_(iy) > by ( |2.6[ ). We also have Imsi^_(iy) > 0, 
so that Im(s? (%)) > 0. This implies Im(l/(3sf + a)) < and so 
indeed by 

-^Resi,_(%) > 0, y>0, 
ay 



which proves part (b) of the lemma. 



□ 



2.5.3 Proof of part (c) 



Proof. Let a < and let V2 be as in Lemma |l.2| (c). Since V2 is a measure 
on iM we have for a; G M, x > 0, 



1 

2 



log(x — z'^)di'2{z), X > 0. 



Hence 



X > 0. 



[x — z^)^ 

Since z/2 < o'2 and the integrand is positive for every z E zM, we have 



dx"^ 



1 



1 



2 

1 d 
2dx 

T d ( Sl 



1 



X — z 



da2{z) 

2d(72{z] 



2 dx \ yfx ) 



where we used the second integral in (2.16). 
Since d'^ = ts\ we see that 



dx"^ 



dx"^ 



Since V3 = 6*2 — ^3, see (2.14), we have for < x < (x*(a))^ 
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0<x< ix*ia)f. 



Since 9i + 02 + 0^ = ^a^, it now also follows that 



(2.19) 



dx"^ 



Recall that 53(2) is the solution of + as = tz with 53(0) = 0. Since 
a < we have that S3 is an odd function which is analytic in a neighborhood 
of 0. Inserting the Taylor series 

00 

= - ^ CfcZ^'^+S \z\ < x*{a), 

with Co = — ^ > into [53(2;)]^ = —as3{z) + TZ and comparing coefficients of 
z, it is easy to show inductively that Ck > for every k. Since O'^(z) = ts3{z) 
with ^3(0) = 0, we then also have 



and 



(2.20) 



d^ 



E 

fc=0 



2A; + 2 



\z\ < X (a) 



2^^3(\/S) = tV fccfcx^' > 0, <x < {x*{a))' 
dx^ ^-^ 



k=l 



since all > 0. The two inequalities (2.19) and (2.20) give the convexity of 
V3{y/x) — If^'^ iVx), which completes the proof of part (c) of Lemma 1.2 □ 



3 Proof of Theorem 11.11 

We basically follow Section 4 of |45j where Theorem 1 1 . 1 1 was proved for the 
case a = 0. However, we need more additional results from potential theory. 



3.1 Results from potential theory 

We use a number of results and notions from logarithmic potential theory. 
The main reference is [78] . Most results in j78j are stated for measures with 
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compact support, while we are also dealing with measures with unbounded 
support, namely the real line or the imaginary axis. Therefore we need 
a number of results from |78) in a slightly stronger form that allows for 
measures with unbounded supports. 

The following theorem is known as the principle of domination, and it 
is stated in \7S\ Theorem II. 3. 2] for the case where // and v have compact 
supports. 

Theorem 3.1. Suppose /x and i' are finite Borel measures with J di' < J dfi. 
Suppose that /x has finite logarithmic energy and that <S'(^) 7^ C. If for some 
constant c the inequality 

(3.1) U^{z) <U''{z) + c 

holds fi-almost everywhere, then it holds for all z G C. 

Proof. Let us first establish Theorem |3 . 1 1 under the assumption that /x has 
compact support, say 5(/i) C Dr = {z \ \z\ < R}. Let i> be the balayage of 
u onto Dji. By the properties of balayage, we then have j du = J dv and 
for certain constant i >0, 

UHz) = U-{z) + i, zGDr, 
UHz) < U''(z) + L z£C. 



Then if (3.1) holds /i-a.e., we find from the equality in (3.2) and the fact 



that S{n) C Dpi that 

U^'iz) <U^{z) + c-£ /x-a.e.. 

Thus by the principle of domination for measures with compact supports, 
see [781 Theorem II. 3. 2], we find < + c — i on C, which in view of the 



inequality in (3.2) leads to If^ < If^ + a on C, as required. 

We next assume that /x is as in the statement of the theorem. As S{fi) 7^ 
C, there is some disk D(zQ,r) = {z £ C \ \z — zq\ < r} with r > that is 
disjoint from S{fi). By translation and dilation invariance of the statement 



in Theorem 3.1 we may assume that D{0, 1) is disjoint from S{fi). We may 
also assume that U'^{0) < 00. 

Let jl be the image of /x under the inversion z 1— t- 1/z. Then fl has 
compact support and straightforward calculations shows that 

(3.3) Uf'iz) = Ui^il/z) - log \z\ [ dfi- U^'iO) 



Two matrix model with quartic potential 



29 



and /(/u) = lifJ-)- If is the image of under the inversion z i— t- 1/z, then 
we hkewise have 



(3.4) 



U^(z) = U^l/z) - \og\z\ / dv - [/''(O). 



Then from (3.1), (3.3) and (3.4), we get 



[/A < [/'^ _log|z| 



dv] + c-C/''(0) + ?7''(0), 



// — a.e. 



Thus C/^ < V^+ci, /i-a.e. where ci = c-[/^(0)+f/''(0) and vi = d/i- 
/ dv)5Q is a finite positive measure with the same total mass as fi. By the 
principle of domination for compactly supported measures ^ and arbitrary 
that we just proved, we find U^^ < U^^ + ci everywhere, which in turn by 
(3.3) and (3.4) leads to < + c. This proves the theorem. □ 



The following result is stated for compactly supported measures in [TSl 
Theorem IV. 4. 5], see also j80] where the result is attributed to de la Vallee 
Poussin [36j. 

Theorem 3.2. Let fj, and v he measures on C with J di' < J djjL, S{fi) ^ C, 
and finite logarithmic potentials and . Suppose that for some c E M, 
we have 



(3.5) 
Let 

Then 



U^{z)<U''{z) + c, zeSifi). 
A = {zeC\ C/^(z) = U''iz)+c}. 



V \a< a* U 

in the sense that v{B^ < /i(-B) for every Borel set B C A. 



Proof. By Theorem |3.1| we obtain from (3.5) that 
(3.6) 



U''{z)<U''{z) + c, 



z G 



It is enough to consider bounded Borel sets B C A. Given such a -B we 
choose R > such that \z\ < R/2 for every z & B. Let fi and z) be the 
balayages of fi and u onto the closed disk Dn := {z £ C \ \z\ < \R\}. By the 
properties of balayage we have, for certain constants ii and £2, 



U^{z) = Uf'iz) + h, U\z) = U'^iz) +i2, ZG Dr. 
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It then follows from ( |3.6| ) that 

(3.7) Ui^iz) <U''iz) + c + ii-i2, z^Dr 



and again by Theorem 
Dr. Equality holds in 



3.1 



the inequality extends to all of C, since S{jl) C 
3.7) for z G Dr n A, and so in particular for z £ B. 
Then by [78, Theorem IV. 4. 5], we have that ^'(-B) < fi-{B). Then also 
^{B) < /u(-B) since B is contained in the interior of Dr and on Dr the 
balayage measures jl and v differ from /i and v only on the boundary BDr = 
{z£C\\z\= R}. □ 



We do not know if the condition S{fi) 7^ C is necessary in Theorems 3.1 



and 3.2 The condition is more than sufficient for the purposes of this paper, 
since we will only be dealing with measures that are supported on either the 
real line or the imaginary axis. 

3.2 Equilibrium problem for 1/3 

Given a measure z/2 < a2 on iM with finite logarithmic energy and J du2 = 
2/3, the equilibrium problem for 1/3 is to minimize 



(3.8) + j {V^{x)-U''^{x))du{x) 

among all measures on M with Jdi^ = 1/3. In case a > we have V3 = 



and then we have that the minimizer 1^3 of (3.8) is equal to 



z^3 = ^z>2 = ^ Bal(zv2,R) 
where 1)2 denotes the balayage of V2 onto M. Then has the density 



dvj _ 1 
dx TT Lm x'^ + \z 



dv2{z) 



and the support of 1^3 is the full real line. This is similar to what is happening 
for the case a = in [45^ Section 4.2]. 

In case a < 0, the external field V3 is positive on the interval {—x* (a) , x* (a)) 
and zero outside. We use this fact to prove the following inequalities for v-^. 



Lemma 3.3. Let he the minimizer for (3.8) among measures on M with 
j du = 1/3. Let c > x*{a), and let 



(3.9) 



Ac = R\{-c,c). 
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Then we have 
(3.10) 
and 
(3.11) 

Proof. The variational conditions associated with the minimization problem 
for (3.8) are 



(^^3) U.> ^Bal(z.2,M) 



(3.12) 



2U''^{x) + V3{x)-U''^{x)>e, xeR, 



where ^ is a constant. Let 1)2 = Bal(z^2)I^) be the balayage of 1^2 onto 
Then U^^ = U"^ on M, so that it follows from (\3.12\) that 

2U'"'{x) = U^^{x) -V3{x) +£, xG 5(1/3)- 

Since V3{x) > 0, we conclude that 

2C/''3(x) < U^^{x) + e, X G 5(1/3), 



By the principle of domination, see Theorem 3.1 (note that the total masses 
of 21^3 and 1)2 are equal), we have that the inequality holds for every x £ C, 



(3.13) 



2U'''{x) < U^Hx)+i, 



X G 



For X G Ac, we have V3{x) = because of the definition (3.9) with 
c > x*(a). Hence 



2^' (x) > U"^ (x) + £ = U^^ (x) + £, 



X e Ar 



because of the inequality in (3.12). Then by inequality (3.13), we find that 
equality holds. Thus 

C {x I 2U'"i{x) = U'^^{x)+l], 



and the inequality (3.10) follows because of (3.13) and Theorem 3.2 



The second inequality (3.11 ) follows in a similar (even simpler) way. Now 
we redefine i>2 as the balayage of U2 onto Ac- 

i>2 = Bal{u2,Ac). 
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Let X £ Ac- Then x G ^(fs) because of (3.10), which has already 



been proved. Then we have, by the property of balayage and (3.12), since 
Vsix) = for X G Ac, 



(x) = U"^ (x) = (x) -e, xeAc = S{u2) 
Thus by another apphcation of Theorem 3.2 we find ( 3.11| ). 



□ 



We note that it follows from (3.10) with c = x*{a), that 
(3.14) (-oo,-x*(a)] U [x*(a),oo) C Sii^s). 



For every c > x*(a), we find from ( 3.11| ) and the explicit expression for the 
balayage onto A^ that 



(3.15) 



dx ~ 



27r a/x^ — 



x^ + Izp 



dv2{z), 



X G 



X > c. 



In the case a < we make use of Lemma 1.2 (c) and Lemma 3.3 



to 

conclude that the support of the minimizer 1/3 is of the desired form. This 
is done in the next lemma. 

Proposition 3.4. Let he a measure on iR with J dv2 = 2/3 and < 

Assume V2 has finite logarithmic energy. Let be the minimizer of (3.8) 



among all measures on M with total mass 1/3. Then the support of is of 
the form 



(3.16) 



S(Z.3)=M\(-C3,C3) 



for some C3 > 0. 

If a < then C3 < x*{a), and i/ C3 > then the density of vanishes 
as a square root at ±03. 



Proof. If a > 0, then S{v'i) = M and we have (3.16) with C3 = 0. 



For a < 0, we use the fact that the external field V3(x) — U^'^{x) is even. 
So it follows from [7H1 Theorem IV.1.10 (f)] that dvi{x) = di>3(x2)/2 defines 
a measure z>3 which is the minimizer of the functional 



I{u) + 2 J (VsiVx) - U'''{Vx))du{x) 

among measures i' on [0, 00) with J du = 1/3. Note that 
(3.17) Sih) = {x^ I X G 5(1/3)}. 
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The external field V3{^/x) — [/'^^(a/x) is convex on [0, {x*{a))'^] by Lemma 
(c), which implies by [73 Theorem IV. 1.10 (b)] that ^(Pgjn [0, {x*{a))'^] 



1.2 



is an interval. 

By (3.14) we already know that x*{a) G S{u3). Then x*(a)^ G S{D3), 
and it follows that 



5(P3)n[0,x*(a) 



[clx*{af] 



for some C3 S [0,3;* (a)]. 



Combining this with (3.14), (3.17), we find (3.16). 

Now suppose C3 = X* (a) > 0. The minimizer 1/3 is characterized by the 
variational conditions 
(3.18) 

f 2U'''{x) = U^'^ix) - V3(x), X G 5(1/3) = (-00, -C3] U [c3,oo), 
I 2?7^3(x) > C/''2(^) _ ^3(3,)^ X G M \ 5(2/3), 



where the inequality on R \ 5(7/3) is indeed strict due to the convexity of 
VziVx) - C/^2(^). Since C3 = x*(q), we have that 2C/^3 = U""^ on 5(2/3) 
and so 2v^ is the balayage of 1/2 onto 5(1/3). Then the density of 1/3 has a 
square root singularity at ±03 and then it easily follows that 



dx 



diesis) 



-00 



as X — )• C3 + . 



This is not compatible with (3.18), since U^'^{x) — Vz{x) is differentiable at 
X = C3. Therefore C3 < x*(a) in case x*{a) > 0. 

Finally, if C3 > then the density of 1/3 vanishes as a square root at ±03 
as a result of the convexity of V3{^/x) — U^^{y/x) again. The proposition is 
proved. □ 



3.3 Equilibrium problem for Vi 

Given 1/2 and 1/3 the equilibrium problem for ui is to minimize 

(3.19) I{v) + j{V{x)-ei{x)-U''^{x))du{x) 

among all probability measures z/ on M. The minimizer exists and its support 
is contained in an interval [— X, X] that is independent of 1^2 ■ 

This can be proved as in |45|, where weighted polynomials were used. 
Here we give a proof using potential theory. 
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Lemma 3.5. Let vi he the minimizer of the weighted energy 



(3.20) + j {V{x)-ei{x))dv{x) 

among probability measures on R, and suppose that 

(3.21) 5(Pi) C 

for some Xi > 0. Let 1/2 be a measure on iM, with finite potential V^'^ and 
suppose that ui is the minimizer for (3.19). Then also 

(3.22) S{ui)c[-Xi,X{\. 

Proof. The variational conditions associated with the equihbrium problems 
for ( [3^ and are 

{2U^^{x) + Vi{x)=ii, xeS{yi), 
\2U^^{x) + Vi{x)>li, xGM, 



(3.23) 
and 
(3.24) 



2U''^{x) + Vi{x)-U''^{x)>li, xGM, 



where li and £1 are certain constants. Combining the relatons (3.23) and 
K2^, we find 



(3.25) 



2U^^{x) - 2[/^i(x) + V^ix) <h-li, X G S(Pi), 
2U''^{x) - 2U'^'{x) - V^ix) <h-ii, X G S{ui). 



Since x 1— )• —U'^'^{x) = Jlog|x — z\dv2{z) is strictly increasing as |x| 
increases, and since S{vi) C [— it follows from the first inequality 
in ( |3.25[ ) that 

2U'^'{x) < 2U''^{x)+ei-£i - V^Xi), X G 5(Pi). 

By the principle of domination, Theorem |3.1[ the inequality holds every- 
where 

2C/^i(x) < 2U'''{x) + £i-£i - U^^Xi), x G C. 

Then for x G 8(1^1) we find by combining this with the second inequality in 
( [335] ) 

-U^^{x) < 2U'''^{x) -2U'''{x) + £i-h < -U^^Xi), 

which implies that |x| < Xi, since x 1— t- — [/^^(x) is even and strictly increas- 
ing as |x| increases. □ 
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3.4 Equilibrium problem for U2 

Given z^i and 1^3 on M with total masses / di^i = 1, / dv^ = 1/3, the equilib- 
rium problem for 1^2 is to minimize 

(3.26) JiU"' +U''^)du 

among all measures on iR with v < (12 and J dv2 = '2./2>. 



Recall that by Lemma 1.2 (b) the density da2{iy) /\dz\ increases as y > 
increases. Then we can use exactly the same arguments as in Lemma 4.5 of 
TS] to conclude that 



(3.27) S{U2) = S{a2) 
and 

(3.28) S{a2 - 1^2) = {-ioo, -102] U [ic2, ioo) 

for some C2 > 0. The proof is based on iterated balayage introduced in 



This proof also shows the following analogue of Lemma 3.3 



Lemma 3.6. Let vi and he measures on M with J dui = 1 and J dv^ 



1/3, and having finite logarithmic energy. Let V2 be the minimizer of ( |3.26 ) 
among all measures v on iM with total mass 2/3 and v < a2- 
Let c > C2 and put 

(3.29) Be = {—ioo, —ic] U [ic, ioo). 
Then we have 

(3.30) (1/2) \b, > \ Bal(i/i + z^3, «M) . 
and 

(3.31) (2^2) |b,< ^Bal(i^i + i/3,Be). 

Proof. This follows from the iterated balayage. Alternatively, it could be 
proved from the variational conditions associated with the minimization 



problems, as we did for Lemma 3.3 We omit details. □ 
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3.5 Uniqueness of the minimizer 



We write the energy functional (1.12) as 



(3.32) E{ui,U2, us) = ^I{ui) + -^IC^^I - 3^^2) + - 2us) 

+ j Vi{x)di^i{x) + j V3ix)du3ix). 

where 2i/i — 3^2 and 1^2 — 2z^3 are signed measures with vanishing integral. 
From this the uniqueness follows as in Section 4.5 of [45j. Note that there 
is a mistake in formula (4.21) of |l5], since the coefficients of I{ui) and 
I{2vi — 3^2) are incorrect. However, the only important issue to establish 
uniqueness is that the coefficients are positive. 

3.6 Existence of the minimizer 

After these preparations we are able to show that the minimizer exists. The 
proof follows along the lines of Section 4.6 of |45j . 

We fix p G (1,5/3). We are going to minimize the energy functional 



(3.32) among all measures vi, U2, 1^3 as before, but with the additional 



restrictions that for certain given X > and K > 0, 
(3.33) S{ui)c[-X,X], 



(3.34) 



du2 

W\ 



< 



K 



for z G iM, \z\ > X, 



(3.35) 



dx ~ xP 



for x £ 



\x\ > X. 



The constants X and K are at our disposal, and later we will choose them 
large enough. 

For any choice of X and K there is a unique vector of measures that 
minimizes the energy functional subject to the usual constraints as well 



as the additional restrictions (3.33), (3.33), (3.33). Indeed, the additional 



restrictions yield that the measures are restricted to a tight sets of measures. 



The energy functional (3.32) is strictly convex, and so there is indeed a 
unique minimizer. 

Our strategy of proof is now to show that for large enough X and K the 
additional restrictions are not effective. 
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Restriction (3.33) This is easy to do for (3.33). Indeed because of Lemma 3.22 



it suffices to clioose 



X > Xi 



where Xi is as in the lemma. Then it is easy to see that (3.33) provides no 
extra restriction. 



Restriction (3.35) The assumption p G (1,5/3) ensures that 



(3.36) 



Cn 



IT 



oo gl-p 



1 + 



ds 



2sin(p7r/2) 



< 1. 



Choose Ep > such that 



(3.37) 

Pick c > such that 

c > x*{a) 
We take X > Xi such that 

(3.38) 

Then X > c and also 
(3.39) 



{l + EpfCp < 1. 



and a2{[—ic,ic\) > 2/3. 



<l + e 



pi 



< 1 + e 



\A > X, 



\x\ > X. 



Assume that f2 is a measure on iM, symmetric around the origin with 



J di'2 and satisfying the restriction (3.34) as well as 1^2 ^ (^2 ■ Let ^3 be the 
minimizer of (3.8) among measures on M with J du = 1/3. Note that we 
do not impose the restriction (3.35). 



Since c > x*{a) we have the inequality (3.15) which due to the symmetry 
of V2 may be written as 



dv-xix) 1 
(3.40) < 



dx 



+ \z\ 



dv2{z) 



G M, |x| > c. 



Let X > X > c. We split the integral in (3.40) into an integral from to 



iX and from iX to ioo, and we estimate using (3.38) 



X 



IX 



y^c2 + |z|2 



c2 Jo + 



dj^2(-z) < 



TT 



max „ , ,„ 

z6[0,iX] X"' + 



I^2([0,iX]) 
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< 



1 + Ep + X2 



and using (3.38) and (3.39) 
1 X 



2 fioo 



IT ^X^ - (? JiX x'^ + \z 



-du2{z) < 



2 , I 12^^2(2) 
JiX X^ + \zY 



X X' 



where the last inequahty holds since V2 satisfies (3.34). This leads to 
1 X 



x"^ — JiX 



x^ + \z\ 



-dv2{z) < 



vr Jo + \z 
tt\x\p Jo 1 + s2 



\dz\ 



ds 



where we made the change of variables z = is, s > 0. Thus by (3.36) we 
have 



/ 9 9 / 2 I I 12 ^^^2^ < 1 |„ , 

7rVx2-c2yix a;2 + |z|2 
In total we get 



X G M, \x\ > X. 



dx 



3tt 



x^ 



xP 



Since (1 + Ep^Cp < 1 and p < 5/3 < 2, it is possible to take now K 
sufficiently large, say K > K^, such that 



37r 



xP 



Then the additional restriction (3.35) is satisfied. 



Restriction (3.34) A similar argument shows the following. Choose vi 



and f3 on M with J dvi = 1, J du^ = 1/3 and satisfying the additional 
restrictions (3.33) and ( |3.35 ). Let 1^2 be the minimizer for (3.26) for 1^2 on 
iR with J dv2 = 2/3 and 1^2 ^<^2- However, we do not impose (3.34). 

Then with the same choice for X > Xi as above (so that (3.38) holds), 
we will find that for K large enough, say K > K2 we have that 



dv2 K 



z eil 



\z\ > X. 
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That is, the restriction (3.34) is satisfied 



This completes the proof of existence of the minimizer for the vector 
equihbrium problem. 



3.7 Proof of Theorem [TT] 

After all this work the proof of Theorem |1.1| is short. 



Proof. Existence and uniqueness of the minimizer is proved in Sections 3.5 
and 3.6 We denote the minimizer by (^i, //2, /^s)- 



The measure /ii is the equilibrium measure in external field Vi — U^'^, 
which is real analytic on M \ {0}. If S'(/ii), then this implies by a result 
of Deift, Kriecherbauer and McLaughlin [32 1 that S{fii) is a finite union of 
intervals with a density that has the form (1.25). If S{a2 — 1^2) then 
Vi — U^'^ is real analytic on M (also at 0) by Lemma 1.2 (b), and again by 



32] we find that S{fj,i) is a finite union of intervals with a density (1.25). 



The conditions (1.23) are the Euler-Lagrange conditions associated with the 



minimization in external field, and they are valid in all cases. This proves 
part (a). 

The statements (1.27) about the supports of /i2 and 0-2 — fj,2 were already 



proved in Section 3.4 see (3.27) and (3.28). The Euler-Lagrange conditions 



(1.26) also follow from this. The fact that p2 vanishes as a square root at 
ibic2 in case C2 > follows as in the proof of Lemma 4.5 of [45j. The other 
statements in part (b) are obvious. 

Part (c) follows from Proposition 3.4, see also (3.18). This completes the 
proof of Theorem |1.1[ □ 



4 A Riemann surface 

The rest of the paper is aimed at the proof of Theorem |1.4[ We assume 
from now on that {V, W, r) is regular, which means in particular that S{ni) 
and S{(J2 — fJ,2) are disjoint. Then by part (a) of Theorem |1.1| we have that 
S{fii) consists of a finite union of disjoint intervals. We use this structure 
as well as that of the supports of the other measures to build a Riemann 
surface in this section. 

We start by collecting consequences of the vector equilibrium problem 
in the form of properties of the (7- functions. These will be used in the 
construction of a meromorphic function on the Riemann surface. 



40 



M. Duits, A.B.J. Kuijlaars, and M.Y. Mo 



4.1 The (^-functions 



The Euler-Lagrange variational conditions (1.23), (1.26), (1.29), can be 
rewritten in terms of the g-functions 

(4.1) gjiz) = j \og{z - s) dfijis), j = 1, 2, 3 

that are defined as foUows. 



Definition 4.1. For j = 1,2,3 we define gj by the formula (4.1) with the 
following choice of branch for log(z — s) with s G S{iJ,j). 

(a) For j = 1,3, we define log{z — s) for s G M with a branch cut along 
(— oo, s] on the real line. 

(b) For j = 2, we define log(z — s) for s £ iM with a branch cut along 
(— oo, 0] U [0, s], which is partly on the real line and partly on the imag- 
inary axis. 

In all cases the definition of\og{z — s) is such that 

log(2: — s) ~ log |z| + i arg z, — 7r<arg2;<7r 

as z ^ oo. 

As a result we have that gi is defined and analytic on C \ (— oo, b^], g2 
on C \ {iR U R-), and 53 on C \ M. 

In what follows we will frequently use the numbers given in the fol- 
lowing definition. 

Definition 4.2. We define 

(4.2) ak = ni{[ak+i, +00)), k = 0,...,N -1, 

and ajy = 0. 

Recall that /ii is supported on UfcLiio^fei^fc] with ai < 61 < 02 < • • • < 
Oat < 67V, and so 

ai = 1 > 02 > • • • > OAf-i > dN = 0. 



The above definitions are such that the following hold. 
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Lemma 4.3. (a) We have 

(4.3) 5i,+ (a;) - 51 = 27ri/ii([x,oo)) for x £ 
In particular 

(4.4) gi^+ - 51 _ = 2Triak, on {hk,ak+i) 



with Ok as in (4.2). Here we have put bo = — oo and oat+i = +oo. 

(b) We have 

92,±{x) = y log |x — s| d;U2(s) ± -Tri, x G M~, 

so that 

A 

(4.5) 

(c) // C2 > then 



92,+ - 92- = -TTZ, on 



(4.6) < 



92,+ {z) - 92 -{z) = -iri - 27rzcr2([0, z]), z G (0, ic2), 
2 

92,+ [z) - 92 -{z) = --TTi + 27ricr2([2,0]), Z £ (-ZC2,0). 



(d) We have 

(4.7) 53,+ (x) -5r3_(x) = 27ri/X3([x,oo)), X e M. 

(e) // cs > 0, i/ien 

53,±(a;) = y log|x - s| dn3{s) ± ^ni, x G (-03,03) 

so that 

(4.8) 53,+ -93- = ^vri, on (-03,03). 



Proof. Parts (a) and (d) are immediate from (4.1). 

Part (b) follows from the definition of §2, the symmetry of the measure 
fi2 on zM, and the fact that / dfi2 = 2/3. 
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For part (c), we note that for z £ 



92,±{z) = -U^^{z) + -7ri±Trifi2{[z,ioo)), 



such that 



(4.9) 92,+ {z) - 92 -{z) = 2mfi2{[z,i(X))) = -iri - 27rz/i2([0, z]) 

since /X2 has total mass 1/3 on iM+. If z G [0,ic2], then ;U2([0,z]) = a2{[0,z]) 



and the first equation in (4.6) follows. The second equation follows in a 
similar way. 

Part (e) follows from the definition of g^, together with the fact that /U3 
is symmetric on M with total mass 1/3, so that //3([c3, 00)) = 1/6. □ 

We have the following jump properties. 

Lemma 4.4. (a) We have 

9i,+ +91- - 92 = Vi -£i, 



(4.10) 



9i,+ + 91 ~ - 92,± = Vi-£i=f -TTi, on S'(/zi) n 



. Re (c/i,+ + gi- - 52) <Vi-£i 
(b) On iM we have 

(4.11) 



on S{ni)nR^ 
on S'(/zi) nR" 
on M\ 5(^1). 



92,+ + 92- = 91 +93, 
Re (52,+ + 92,-) > Re {91 + 93) 



on S{a2 - ^2), 

on iM \ 5((T2 — fJ-2)- 



(c) We have on 



(4.12) 



93,+ + 93,- -92 = V3, 



93,+ + 93,- - 92,+ =V3^ -TTi, onS{fi3)n 



. Re (53,+ + 93,- - 92) <V3, 



on ^(/ia) nM^ 
on 5(^3) nM" 
on ]R\ 5(/i3). 



Proof. Part (a) follows from (1.23). Indeed, we have 

2C/''i(x) = -5i_+(x) - gi _(x), for X G 
\-92{x), 



for X G 
~92,+ {x) + l^ri, for x G 
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so that ( |1.23D indeed leads to ( |4.10D . 
For part (b) we note that 

Also because of symmetry of gi and 33 around 0, 



f/M3 



for z G iM=^. 



Using this in (1.26) we obtain (4.11). 

Part (c) follows from (1.29) in the same way that we obtained part (a) 
from ([T^. □ 

In what follows we also use the derivatives of the ^(-function, which we 
denote by Fi, F2, F^. 



i = 1,2,3, 



Definition 4.5. We define 

(4.13) F,{z) = g'jiz) -- 

which is defined and analytic for z £ C \ S{fij). 

The jump of Fj gives us the density of fij, since we have 



dfij{s) 



z — s 



(4.14) 



dm 

dx 
dn2 
dz 
dH3 
dx 



1 

'2TTi 

1 

'2TTi 

1 

'2^ 



(^1,H 


_(x) 


-^1,- 


-(^)), 


X G I 




_(z) 


-F2,- 


(^)), 


z G i 


(^3H 


-(x) 




-(^)), 


X G I 



4.2 Riemann surface TZ and ^-functions 

We construct a four sheeted Riemann surface IZ in the following way. Four 
sheets IZj defined as 



(4.15) 



(ni = c\s{m), 

7^2 = C\(5(/il)U5(a2-^2)), 
7l3 = C\(5(a2-/i2)U5(^3)), 

l7^4 = c\5(A^3), 
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are connected as follows: TZi is connected to 7^2 via S{fj,i), IZ2 is connected 
to 7^.3 via S'((T2 — /U2) and IZ^ is connected to IZi via S{^^). Every connection 
is in the usual crosswise manner. We compactify the Riemann surface by 
adding a point at infinity to the first sheet IZi, and a second point at infinity 
which is common to the other three sheets. 

The genus of the Riemann surface is (recall that S{^i) consists of N 
intervals and recall the classification of the cases in Section |1.5|) 



(4.16) 



1, 



in Cases I, II, and III, 
in Cases IV and V. 



Using the functions Fj defined in Definition 4.5 we define the ^ functions. 
Definition 4.6. We define functions ^j, j = 1,2,3,4 by 

(^^ = V'-Fi, onTZi, 

C2 = Fi-F2 + 0[, on 7^2\«, 

= F2 - F3 + onTZ3\(R\JiR), 

U4 = i^3 + f3, on 7^4\M. 



(4.17) 



We first prove that ^2, '^3, and ^4 are analytic on their full respective 
sheets. 



Lemma 4.7. (a) We have ^2, 
analytic extension to 7^2- 



on {—ic2,ic2), and so ^2 has an 



(b) We have ^3^-(_ = ,^3^_ on {—ic2,ic2) and on (—03,03), and so ^3 has an 
analytic extension to TZ^. 



(c) We have ^4^. 
to 7^4. 



^4^_ on (—03,03), and so .^4 has an analytic extension 



Proof. On {—ic2,ic2) we have by (4.6) and (4.13) that 

d 



F2,+ {z)-F2, 



dz 



{g2,+ {z)-g2-{z)) 



-2m- 



.da2 

dz 



= r(si,+ (z) - si _(z)), z € {-ic2,ic2), 
see (2.7). Since rsi = 9[, we get 

i^2,+ - 6*1,+ = F2- - 0[_, on (-io2,io2). 
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and also 



6 



on {-ic2,ic2) 



since -Fi is analytic on the imaginary axis. This proves part (a). 
We also have 

si^+{z) - si -(z) = S2-{z) - S2,+ {z), z G iM. 

Then above argument also shows that 

C3,+ = 6-, on(-ic2,ic2) 

since -Pi and F3 are analytic on the imaginary axis. 

If C3 > (which can only happen if a < 0), then P3 is analytic across 
(—03,03). Both ^2 and ^3 are analytic across (— x*(a), x*(a)), and so a 
fortiori across (—03,03), since 03 < x*{a). This proves part (c) and the 
remaining statement of part (b). □ 

We continue to denote the analytic extension by ^j, j = 2, 3, 4. 

Proposition 4.8. The function 

4 

^ U 7^,■ ^ c 

i=i 

given by ^(z) = (,j{z) for z G TZj extends to a meromorphic function (also 
denoted hy ^) on TZ. The meromorphic function has a pole of order deg F — 1 
at infinity on the first sheet, and a simple pole at the other point at infinity. 

Proof From ( |4.10[ ) and ( |4.13[ ) it follows that 

(4.18) Fi,+ (x) + Fi,_(x)-F2(x) = Fi'(a;), x G S{^ll), 

since F2 is analytic on M \ {0}, so that F2,± = ^2- Since Vi = V — 61, we 
obtain by the definition (4.17) that 

(4.19) Ci,±{x)=C2,t{x), x€S{f,i). 



From the first equation in (2.10) and the fact that 02 > y*{a), we obtain 
that 

6'i,± = 6*2,^5 on 5((T2 — ^2) = (— ioo, —^02] U [zo2,oo). 
Then we obtain from ( 4.11[ ) that 

(4.20) F2,+ iz)+F2,-{z) = Fi{z) + F3{z), 2 G 5((72 - ^^2), 
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and it follows from (4.17) that 

(4.21) 6,±(^)=6,t(^)> zeS{a2-fi2)- 
From (4.12) we similarly find 

(4.22) F3,+ (x) + F3,_(x) - F2(x) = V:^{x), x E 5(^3). 
We next claim that 

(4.23) V3{x) = e2,±{x)-d3^^{x), xeR. 



Indeed, for 2; G {-x*{a),x*{a)) we have by ( |2.15[ ) that V^ix) = 62{x)-6^{x) 
and (4.23) holds, while for 2; G M with |2;| > x*, we have by (2.10) and (2.15) 
that both sides of (4.23) are equal to zero. Using (4.23) in (4.22) we obtain 
from the definition ( 4.17[ ) that 

(4.24) 6,±(^) =^4,t(^), XG5(;U3). 

The analyticity of across the three cuts S{iii), S{a2 — fi2) and S{fj,3) 
is now estabhshed by ( |4.19D, (|4.21D , an d (|4.24D . 
As z — >• 00, we have by (4.17) and ( 4.13[ ) that 

Ci{z) = V'{z) + 0{z-'), 

which implies that ^ has a pole of order deg V — 1 at the point at infinity 
on the first sheet. 

From ( |4.17[ ), ( |4.13[ ) and Lemma [2^4] it also follows that for j = 2,3, 4, 

e,(z)=a;V/V/3 + 0(z-i/3)^ 



as z 



00, 



where the value of k depends on j and on the quadrant in which 2; — )■ 00. 
Since the other point at infinity (which is common to the second, third 
and fourth sheets) is a double branch point, we have that z~^^^ is a local 
coordinate, so that ^ indeed has a simple pole. □ 



Remark 4.9. It follows from Proposition 4.8 that ^j, j 
solutions of a quartic equation 



1,2,3,4 are 



(4.25) 



+ P3{z)f + P2iz)f + piiz)^ + poiz) = 



with coefficients that are polynomial in z. This algebraic equation is known 
as the spectral curve [13]. The degrees of the polynomial coefficients are 
determined by the degree of V. 
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Using this fact, we indicate how to remove the condition that S{fii) 
and 5 ((72 — ^2) are disjoint in part (a) of Theorem This condition 

was included in order to be able to conclude that 5'(/Ui) is a finite union of 
intervals. In case the condition does not hold, we can now argue as follows. 

Given e > we modify the vector equilibrium problem by requiring that 
^(/ii) C (—00, — e] U [e,oo). Then the proof of existence and uniqueness of 
the minimizer follows in the same way as in Section [3} Let us denote the 
minimizer by (/xf , //Ij A'l)- Because the external field Vi — ig real analytic 
on M \ (— e, e) we have that the support of /xf is a finite union of intervals. 
The further structure of the minimizers is the same, which means that we 
can construct a Riemann surface with a globally meromorphic function on 
it, in the same way as we did in this section. The only difference is that 
the meromorphic function may have a pole in ite. The algebraic equation 



(4.25) has coefficients that are rational in z, with ±e as the only (simple) 



poles, and so we may write the spectral curve as 

(4.26) (Z^ - £2)^4 ^ ^^(^)^3 ^ ^^(^)^2 ^ . 







with polynomial coefficients qj{z), whose degrees is determined by the degree 
of V. In particular, the degrees do not depend on e. 

As e — )• we have that fij converges weakly to fij for j = 1,2,3, and 
it is not difficult to show that the spectral curve (4.26) has a limit as well. 
Then V-l{z) — J '^^^^^^ is the solution of an algebraic equation and therefore 
the support of m is a finite union of intervals. 



4.3 Pro 



(perties of the ^ functions 

.17) we find that 



From (4.14) and the definition (4. 

(a, 



dm 


1 


dx 


~ 27ri 


dm 


1 


dz 






1 



X £ 



[x) - (.i-ix)) , 

;^)-6,-(^))-^ K+(^)-^'i,-(^)) 

t? = ^(e..(x)-6.(x))-^(^^,,(x)-.^^_(-)) 

^ - (x)-e4,-(^))+2^(^3,+ (^)-^3,-(^)), 



z e il 



X G 



Since ^ is an odd function with a simple pole at infinity on sheets 2, 3 
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and 4, there is an expansion in the local coordinate z 



C^lZ^^^ + CiZ ^^^ + C3Z ^ + C5Z 5/3 + .. 



with certain real constants cj. Keeping track of the principal branches, we 
obtain the following for the asymptotic behavior of ?3 and ^4 as z — t- 00. 
Note that the structure of the formulas in terms of the factors w and uj'^ is 
the same as in Lemma EiS] and Lemma [231 



Lemma 4.10. We have as z —t- 00 





/■ 

c_ 






C5Z-5/3 + 




in 


luiv, 


= < 


C- 


-lUZ^/^ + ClW^Z 




'1 + CsWZ" 


5/3 + ... 


in 


II, 




c_ 

V 


.itJ^Z-'-/^ + ClOJZ 


-1/3 + C32;- 


"1 + C^Uj'^Z 


-5/3 _^ . . . 


, in 


III, 




/■ 

C- 


^lOJZ^/^ + ClW^Z 


-1/3 + C3Z- 


'1 + cswz- 


5/3 + ... 


in 


I, 


= < 






+ C3Z"1 + 


C52;~5/3 + 


* * * ) 


in 


II u ///, 






_itj^z-'-/3 + ClWZ 


-1/3 + C3Z- 


'1 + csw^z 


-5/3 + . . . 


, in 


IV, 






lOj'^Z^^^ + ClWZ" 


-1/3 + C3Z- 


1 + C^Uj'^Z 


-5/3 + . . . 


, in 


lull. 


I- 




-1/3 + C3Z~ 


1 + C5a;z~ 


5/3 + ... ^ 


in 


III U IV. 



with constants 
(4.28) c_i 



r4/3 



Cl 



.« 2/3 

3 ' 



C3 = 1/3. 



Proof. The form of the asymptotics follows from the definition (4.17) snd 



Lemmas 2.3 and 2.4 Since /xi is a probability measure with compact sup- 
port, and symmetric with respect to the origin, we have 



(4.29) 
From 
(4.30) 



1 



Fi{z) = - + 0{z-'). 



F2{z) = e[{z) - Uz) + Fi{z) = Tsiiz) - 6(^) + Fi{z) 



we see that F2 has a series representation around z = 00 in powers of zi/3. It 
should start with F2{z) = i; 



+ • • • since ^2 has total mass 2/3, see (4.13). 



Then the coefficients for 2:1/3 and z i/3 should vanish which by Lemma 



2.3 



and (4.30) leads to the expressions in (4.28) for c_i and ci. The coefficient 
of z 1 should be 2/3, which by (|4.29D and (4.30|) leads to C3 = 1/3. □ 
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Combining (4.29) and (4.30) with Lemmas 2.3 and 4.10 we also find that 



there is a real constant 
(4.31) 



such that 



(4.32) 



r2^-l_Cz-5/3 + 0(z-2), 



|z"i-Ca;z-5/3 + 0(z-2), 



ZGIUIV, 
z G II, 
z G /// 



as z — >• oo. Then by (4.17), and Lemmas 2.4 and 4.10, we also get 
(4.33) F3(z) 



1^-1 



Z- - + Coj^z-^/^ + 0{z~^), ZGIUII, 



1 ^-1 



+ Ca;z-5/3 + 0(z-2), 



zelllu IV, 



as z — 7- oc, with the same constant C. Using this and Lemma 2.4 in (4.27) 
we find that 



(4.34) 



^ = ^6-1^1-5/3 + 0(^-2), as\z\^oo,zeil 
\dz\ 2tt 

^ = ^Cx-5/3 + 0(x-2) 

dx 2tt ^ ' 



), as X — oo, X € M. 

which gives the precise rate of decay of the densities of /i2 and along the 
imaginary and real axis, respectively. It also follows from (4.34) that the 
constant is positive, C > 0. 

Furthermore, after integration, we find from (4.29), (4.32), and (4.33) 



giiz) = \ogz + Oiz-') 

logz+ |Cz-2/3 + o(z-i), 
g2{z) = { f logz + |Ca;z-2/3 + 0(z-i), 



3 ^& ^ , 2 ^ 

|logz + |Ca;2z-2/3 + o(2-i), 

llogz-lCuj'^z-'^/^ + 0{z-^), 
llog z - ICuz-^/'"^ + 0(z'^), 



ZGIUIV, 
z G //, 
z G ///, 

z e lull, 
zelllu IV 



as z ^ oo. It may be verified from (4.1) that the constant of integration 
indeed vanishes. 
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4.4 The A functions 

The A functions are primitive functions of the ^ functions (that is, Abehan 
integrals). They are defined as fohows. 

Definition 4.11. We define 

( Xi = V-gi-h, 

h = 91 - 92 + Oi, 

h = 92 - 93 + 6*2, 

I A4 = C/3 + O3. 



(4.35) 



It is convenient to consider each Xj function as defined on TZj with an 
extra cut to ensure single- valuedness. Thus, Ai is defined and analytic on 
C\(-oo,6jv], A2 on C\((-oo,6Ar]UiM), A3 onC\(MUiM), and A4 onC\M. 

The A- functions satisfy the following jump conditions that follow by 



combining the jumps for the ^-functions given in Lemmas 4.3 and |4.4| with 
the properties (2.10) of the functions 9j. 

Lemma 4.12. 




0, 



A2,T = ^FgVri, 



-2Triak, 

4 

2'iriak - -vri, 



on 5(;Ui) nM+, 
on 5(^1) nR". 

on (6fc,afc+i) n : 



on ibk,ak+i)n : 



where bo 
0,...,iV. 



-00, oat+i = +00, and ak is given by (4.2) for k 



(b) On the imaginary axis we have 
(4.39) A2,± - As,^ = 

and 

f 

A2,+ — A2,- 



on S{a2 - II2) 



(4.40) 



A 



3,H 



A3 



^3 ' 
,2 . 

3 ' 



on {—ic2, ic2) n ' 
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(c) Finally, 
(4.41) 
and 
(4.42) 



\3,± 



0, 



on (c3,oo), 
on (-00, -C3), 



A3, 
A4, 



on (-C3,C3). 



Proof. The equations follow by combining the jumps for the (7- functions 
given in Lemmas 4.3 and 4.4 with the properties (2.10) of the functions Oj. 
Only (4.40) requires more explanation. 

get for z e 

\TTi + 27rm2([0, z]) + (ei,+ (z) - ^l,-(^)) 



From (4.6) and (4.35) we get for z G (0,ic2) 

2 



A2,+ (2;) - A2 -(2;) 



where we used the fact that gi and (73 are analytic on the imaginary axis. 
Since (T2 has density (2.7) and 9[ = rsi, we have 

da2 



dz 



{z)dz 



1 

1 

27ri 



9i,+ (z)-^i,_(O)+0i,+ (O)) 



and the first equality in (4.40) follows for z G (0, zc2) since 61 -{Q) = 9i^j^{Q). 
The other equalities contained in (4.40) follow in the same way. □ 

The A functions also satisfy a number of inequalities. 

Lemma 4.13. We have the following inequalities 



(4.43) 
(4.44) 
(4.45) 



Re(A2,+ - Ai _) < onM.\S{iii), 
Re(A2,+ - A3,_) < on iM\5(cJ2 -^2), 
Re(A4,+ - A3-) < onM.\S{iiz). 



Proof. The inequalities follow from the inequalities in Lemma 4.4 



For example, to obtain (4.43) we note that by (4.35) and (2.14) 
A2,+ - Ai,_ = (gi - 52 + ^1)+ -{V-gi- h)- 
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= gi,+ + 91- - 92,+ -Vi + h 
which indeed has a negative real part on M \ S{p,i) because of the inequahty 



in (4.10). The inequahty is strict because of the regularity assumption. 



The other inequalities (4.44) and (4.45) follow in a similar way 



□ 



The asymptotic behavior of the A-functions follows by combining the 



definition (4.35), which we state for ease of future reference. 



Lemma 4.14. We have as z —)■ oo 



(4.46) 


Ai(z) 


= V{z) — logz — 


il 


+ 0{z-^), 










^llo, 


[Z 


-lCz-^l^ + 0{z-^), 


z G I u /y, 


(4.47) 


A2(^) 


= 0i{z)+< 




[Z 


-iCa;z-2/3 + 0(z-i), 


z G //, 








[llog 


[Z 


-|Cl^2^-2/3 + 0(z-1), 


z G ///, 










[Z 


-iCa;z-2/3 + 0(z-i), 


z G /, 


(4.48) 


A3(^) 


= e,iz) + < 




[Z 


-|Cz-2/3 + 0(z-1), 


z£ IIU III, 










[Z 


-|Cl^2^-2/3 + 0(z-1), 


z G IV, 


(4.49) 




= e^iz) + < 


filog 
1 5 log 


z 
z 


-|Cc^2^-2/3 + 0(z-1), 
- |CtJZ-2/3 + O(z-l), 


z e lull, 

zelllU IV, 



where C is the constant from (4.31). 



5 Pearcey integrals and the first transformation of 
the RH problem 



Now we come to the first transformation of the RH problem (1.36) which as 
in [l5] will be done with the help of Pearcey integrals. The present setup is 
however slightly different from the one in [13], since we will construct a RH 
problem for the Pearcey integrals that is n-dependent. 

5.1 Definitions 



We note that the weights from (1.35) can be written as 



(5.1) 



e-"^("-Vo,n(x) 



wi,n{x) = (nr)-ie-"^(^Vo,n(^), 
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where 

(5.2) Po,n{x) 



-n(W(s)-TXs)^^ 



W{s) 



1 4 « 2 

4^ +r 



is a Pearcey integral that satisfies the third order linear ODE 

/r o\ III . 2 2 I 34 n 

(5.oj "p + n T ap — n T xp = U. 

Other solutions to the same ODE are given by similar integrals. 
Definition 5.1. For j = 0, . . . , 5 and n gN, we define 



(5.4) 



-n{W{s)-TZs)^^ 



zeC, 



where the contours Fj are 



(5.5) 



To = (-00,00), 
r2 = (ioo,0] U [0,-00), 
r4 = (-ioo,0] U [0,00), 



Ti = (iioo,0] U [0,00), 
Ts = (-ioo,0] U [0,-00), 



or homotopic deformations such as the ones shown in Figure [XH and with 



the orientation as also shown in Figure 5.1 



All Pearcey integrals (5.4) are entire functions in the complex plane. 
Certain combinations of these functions are used to build a 3 x 3 matrix 
valued P„ : C \ (M U zM) C^^^ as follows. 

Definition 5.2. In each of the four quadrants (denoted I, II, III and IV) 
we define 



(5.6) 



^PO,n 
P'o,n 
\pln 
^PO,n 
P0,n 
\P0,n 
/P0,n 
P0,n 

/P0,n 
P0,n 
\p'in 



-P2,n 
-P2.n 
-P'in 
-Pl,n 



-P5,n 
~P'h,n 
~Ph,nJ 
-P5,n^ 



-Pl,n 
'Pl,n 
-PA,n 
~PA,n 
-Pin 
-P3,n 



P5,n 
-Pln^ 
-Pb,n 

-pin 

-Pi 
-P5,1 



-P3.n -P5,: 



in I, 



in II, 



in HI, 



in IV. 



-Pin 



pin/ 
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^2 1 

^ I 


5 














] 




v\ 





Figure 5.1: Contours Fj in the definition of tfie Pearcey integrals 



Tfien is indeed defined and analytic in C \ (M U iM) with the follow- 
ing jump properties on the real line (oriented from left to right) and the 
imaginary axis (oriented from bottom to top). As usual the orientation de- 
termines the + and — side of the curve, where the + side is on the left and 
the — side on the right. 



Lemma 5.3. T/ie matrix valued function Pn satisfies the jump relations 
(5.7) 



'1 0^ 

Pn,+ {z) = Pn,-{z) \ 1 
.0 -1 1, 



'1 -1 0> 

(5.8) P„,+ (z) = P„,_(z) I 1 

.0 1, 



z G 



z e 



Proof. The jumps easily follow from the definitions. It follows for exam- 
ple from the definition of the Pearcey integrals and the orientation of the 
contours Tj , see Figure 15.11 that 



-P2,n = -P3,n + P5,n- 



This relation implies the jump relation (5.7) for z G M"*" in view of the 



definition of P„ in the first and fourth quadrant given in (5.6) 
The other jumps are proved in a similar way. 



□ 
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5.2 Large z asymptotics 

The large z behavior of Pn{z) is obtained from a classical saddle point anal- 
ysis of the integrals (5.4) that defines each of its entries. Recall that the 
saddle point equation (2.2) has the three solutions si{z), S2{z) and S3{z), 
The value of —W{s) +tzs at the saddle Sj{z) is denoted by 



2.1 



see Section 



9j{z), as in Section 



2.2 



Lemma 5.4. We have as z ^ oo, 



(5.9) Pn{z) 



/27rn 



' n 



^-4/3^-1/3 







-2/3 







-1 +0{z 



-4/3x 







X < 











— UJ \ 


(i 




-1 




-1 


■] 






-UJ 




—u 




/- 






1 




'UJ 




-1 




1 




-1 


v- 


-LO 




1 




UJ^ 


/- 


-OJ 




1 


u: 






-1 




1 


1 


1 


V- 


OJ^ 




1 


UJ 1 


/l 




-UJ 




uj-\ 




1 




-1 






I- 


\1 








I) 





pne(2) 



^ne{z) 



z in /, 



Z in II, 



z in III, 



z in IV, 



where G is defined as in (2.11). 

Proof. The proof is a tedious saddle point analysis for all integrals that 
define Pn in the respective quadrants, see also [21j. The ODE (5.3) can be 
used to find the form of the asymptotic expansion. Indeed, putting p = e 
in (5.3) we obtain the following differential equation for f = F' , 

f 



nF 



r^z + r'af+^ff + ^f" = 0. 
This nonlinear ODE has solutions / with expansions 
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fc^4/3^1/3 



CO T 



+ 



a^2fe^2/3^-l/3 

3 



3n 



a 



\81 9n 



k -2/3 -5/3 



UJ T 



+ 0{z 



-7/3N 



k = 0,1,2, 



as z — >• oo in one of the quadrants. Here uj = e^'^*/^, as before. Thus after 
integration and by Proposition 2.4, we have for a certain j = 1, 2, 3, 



F{z) 



1 

3n 



1 



log z + - log C + -^w'=(rz)-2/3 + 0(^-4/3) 



n 



6n 



where C is a constant, and so there are solutions of the Pearcey ODE ( |5.3[ ) 
that behave like 

(5.10) p{z) = ^ «^fc(^^)-2/3 ^ 0(z-4/3)) 

as z — )• oo in one of the quadrants. 

Each of the functions po^n, ■ ■ ■ ,P5,n that appears in the definition of Pn 
in a certain quadrant has an asymptotic expansion as in (5.10). We have to 
associate with each such function a value of j and the corresponding value 
of k. To do this, we have to perform a saddle point analysis on the integral 
representation (5.4). 

The saddle point equation 

(5.11) W'{s)-Tz = s^ + as-Tz = 



for (5.4) has three solutions Sj{z). It turns out that in each quadrant and 
for each j = 1, 2, 3, we have that Sj{z) is the relevant saddle for the Pearcey 
integral that defines the function in the jth column. Thus the initial contour 
can be deformed into the steepest descent path through Sj{z), or into a union 
of steepest descent paths with Sj{z) as the determining saddle. 

Then by a classical saddle point analysis, see e.g. [73], we obtain the 
following asymptotic behavior 



p{z) 



2tt 



(5.12) 



nW"{sj{z)) 
nW"{sj{z))' 



-n{W{sj {z))—TzSj(z)) 



(l + 0{z- 



2/3- 



^nOj (z) 



l + 0(z-2/3; 



as z — )• oo. Since W"{s) = 3s^ + a, and Sj{z) = 0{z^^^) by Lemma 
3 behavior 

p{z) = 



find from (5.12) the behavior 



2.3 



we 




-lne,(z) 



{l + 0{z 



2/3^ 
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which then by (5.10) takes the form 



V 3n V 6 / 

for a certain value of k, depending on j and depending on the particular 
quadrant. In the first quadrant, for example, we have fc = j — 1 as follows 
from Lemma |2.3[ 

We may differentiate the expansions and we obtain for the derivative 

V on V 6 / 

= ± l^nre'^eA^^ ( 1 - -u;'^{tz)-^/' + 0{z-^/^)) 

V 3n V 6 / 

and for the second derivative 

Y Stz- V 2 / 

In each quadrant we have to take the correct sign (determined by the 
orientation of the contours) and the correct value of k. The formulas can 
then be written in the form (5.9). □ 

Let us define the constant matrices Aj as follows. 
Definition 5.5. We define 



(5.13) Ai = ^il -1 -1 , A2 = ^\ -I -1 -1 




(5.14) ^3^^ -1 -1 1 , -1 1 




1 -u^' 
1 -1 

1 -UJ 



'IjJ^ CO 



U) up' 



The prefactor ^ is chosen such that det Aj = 1 for j = 1, 2, 3, 4. Then 
we can reformulate Lemma 15.41 as follows. 



Corollary 5.6. We have as z —)• oo in the jth quadrant, 

(5.15) P„(z) = P„,o(/ + 0(^"'/'))diag(z-i/3 i ^1/3) AT*e"®(^) 

where Pnfi is the invertible matrix 

'n-^T-'^l'^ 

(5.16) P„o = \/2™ri( 10 

-fnr2/3 nr^/^ 
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The asymptotic formula (5.15) will be the most convenient to work with 
in what follows. Note that we still have an error term / + 0{z~'^^^) which 
is somewhat remarkable. 



5.3 First transformation: Y X 



With the 3x3 matrix- valued Pn we can now perform the first transformation 
of the Riemann-Hilbert problem. Recall that Y is the solution of the RH 
problem (1.36). 



Definition 5.7. We define X by 
(5.17) ^W=(J 



1 

L>„P^*(z)e"®(^) 



for z in the jth quadrant. Here Cn and Dn are the constant matrices 
(5.18) Cn = PloD-\ Dn = diag{l ut {nrf), 



withPno given by (5.16), and Pn is given by i5.6h, andO is given by (2.11). 



The matrices in the right-hand side of (5.17) are 4x4 matrices written in 
block form, where the right lower block has size 3x3. The transformation 
(5.17) does not affect the (1,1) entry. The factor e"'^{z) is included in 
the definition of X(z) in order to simplify the asymptotic behavior of X. 
However, it will complicate the jump matrices, as we will see. 

The matrix valued function X is defined and analytic in each of the four 
quadrants. 

The asymptotic behavior for X is as follows. 



Lemma 5.8. We have 



(5.19) X(z) = (/ + 0(z-2/3)) 



n I 1 

Z 3^3 










Z 3 





\ 



_n_l 
Z 3 3 



1 

A, 



as z ^ oo in the j th quadrant with matrices where the matrices Aj are given 
by dslsl) and (Isll). 



As in ^45j we have that the asymptotic formula ( 5.19 ) for X is not uniform 
up to the axis. 
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Proof. By Corollary 5.6 and the definitions (5.18) of C„ and Dn we have 

\ 

aZ?nP-*(^)e"®(^) = + 1 

^ Moo Z-V3/ 



as z — )• oo in the jth quadrant. Using this in (5.17), together with the 
asymptotics of Y as given in (1.36), we obtain the lemma. □ 

The jumps for X on the real and imaginary axis are as follows. 

Lemma 5.9. We have 

X+{z) = X^{z)Jxiz), z£Rum, 

where the jump matrices Jx are given as follows. 

• On the real line we have for x G (— oo, —x*{a)] U [x*{a),oo) 

(5.20) 



Jxix) 



\ 



1 

e"(^2,+ (2.')-e3,+ (a;)) \ 

yO Qn{e-i^j^{x)-Qi^^{x)) J 



and for X G (— x*(a)) (only relevant in case a < Oj, 

/I (,-'^{V{x)-ei{x)) Q \ 

10 

1 e~"(^2(a;)-e'3(3;)) 

Vo 1 / 



(5.21) .Jx{x) 







On the imaginary axis, we have for z G {—ico,—iy*{a)\U[iy*{a),ioo), 

(5.22) Jx{z) 

and for z G (—iy* (a) , iy* (a)) (only relevant in case a > 0), 

(5.23) Jxiz) 



/I 

eHSi.+iz)~e2.+(z)) Q 

Vo 1/ 



/I 





0\ 

1 

Vo 1/ 
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Proof. The jump matrix for x G M is by (5.17) and the jump condition in 

(pel) 



Jxix) 



1 

e-"®-(^)P^_(x)D-i 



(5.24) 



fl Wo^n{x) Wi^nix) W2,n{x)\ 
1 

1 

yo 1 / 

1 

1 {wo,n{x) wiA^) t^2,n(x))l)„P„-^(a;)e"Q+(-)\ 

e-"®-(^)P*^_(x)P^^(x)e"®+(^) / 



The row vector in the right upper corner of (5.24) is by (5.1) and (5.18) 

(5.25) {wo,n{x) wi,n{x) t/;2,n(x)) i?„P„-^(x)e"®+(") 

= e-'^^(^) (po,n(x) p'o,„(x) pI,{x)) p-^(x)e"®+(-) 



Since {po,n{x) Po„(x) Po„(2;)) is the first column of Pnix), see (5.6), it 
follows that ( [5^ is equal to e-"^(^) (l O) e"®+(^) which by ( pUj ) leads 
to 

(5.26) {wo,n{x) wi,n{x) «;2,„(x)) i?„p-^(x)e"®+(") 

^ ^g-„(v{x)-ei(x)) Q o) . 



This leads to the first row of the jump matrices (^5.20)-(5.21 ). 



To evaluate the 3x3 block in the right lower corner of (5.24) we note 



that we have by (5.7) 



Pi-{x)P~X{x) = (P„-i(x)P„,+ (x))" 



1 0> 
1 1 
.0 1, 



so that 



'1 0^ 
1 1 1 e"®+(^') 



(5.27) e-"®-(^')p*^_(x)P„-^(x)e"®+(^') = e-"®-^^) 

\0 1^ 

Then using (2.12) to write e-"0-(^) in terms of e-"0+(^)^ and then using 
the explicit expressions (2.11) for 0, we see that (5.27) indeed reduces to 
the 3 X 3 right lower block in ( |5.20[ )-( [5^ . 

The proof of the other expressions for Jx follows in a similar way. □ 
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5.4 RH problem for X 

To summarize, we have found the fohowing RH problem for X 
(5.28) 

' X is analytic in C \ (M U iR), 



X+ = X^Jx, 



on 

2 , 



X{z) = {I + 0{z- 3)) diag 



n — 1 n 

Z 3 2 3 



n + 1 
Z 3 



1 

A, 



as z — >• oo in the jth quadrant, 



where Jx is given by ( 5.20 )-( 5.23). 



Each of the jump matrices Jx is nontrivial only in certain 2x2 blocks. 
The nontrivial blocks are triangular and assume one of the forms 



1 * 
1 

with a real off-diagonal entry *, or 



or 



or 



1 

* 1 



* 
1 * 



with oscillatory diagonal entries of absolute value 1 . The first form indicates 
that an external field is acting and the second form indicates the presence 
of an upper constraint. In this way we can already see the connection with 
the vector equilibrium. 

Let us examine this in more detail. 



Jump Jx on the real line The jump matrix Jx on the real line, see 
(5.20) and (5.21), takes the block form 

where {Jx)i and (Jx)3 are 2x2 matrices. 
We have 



(5.29) 
where 



I ^-nVi{x) 
1 



X G 



Vi{x) = V{x)-9i{x) 
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is indeed the external field that acts on the first measure in the vector 



equilibrium problem, see (2.14) 
)re, we ha 



Furthermore, we have by (|5.21|) 
(5.30) 
with 



1 



X G (— x*(a), x*(q)), 
V3{x) = 02{x) - 03{x) for X G (-x*(a),x*(a)), 



which by (2.15) is indeed the non-zero part of the external field V3 that acts 
on the third measure in the vector equilibrium problem. The external field 
V3 plays a role only in case x*{a) > 0, that is, in case a < 0. 

The right lower block in (5.20) has oscillatory diagonal entries. We define 

Tps by 



(5.31) 
so that 

(5.32) {Jx)3{x) 



,ni/'3,+(a:) 





1 



X G M, |x| > x*{a). 



Then 'ip3^± is purely imaginary for |x| > x*{a) with 

^V'3,±(x) = ^[e2,+ (x)-e2,-(^)] 

= r(s2,+ (x) - S2-{x)) = 2iTlms2,+ ix), 

which is purely imaginary with positive imaginary part. Thus we can asso- 
ciate with ■03 a measure CJ3 on M by putting 



(5.33) 



das 
dx 



1 d 
27ri dx 



{62,+ {x) - 6'3_+(x)) = - Im S2,+ {x), for x G 



vr 



Then 5(^3) = M if a > and Sia^) = R\{-x*{a),x*{a)) if a < 0. 

Because of the upper triangular form of ( 5.32[ ) it will turn out that (T3 
acts as a lower constraint on the third measure in the sense that 



(5.34) 



Ms + 0-3 > 



and then we could allow signed measures /i3 in the vector equilibrium prob- 
lem. However, in this more general vector equilibrium problem we would 
still find //3 > so that the constraint (5.34) does not play a role after all. 
We will not use the measure CJ3 anymore. 
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Jump Jx on the imaginary axis The jump matrix Jx on the imaginary 
axis, see (5.22) and (5.23), takes the block form 



'1 0^ 
Jx= {0 {Jx)2 I , 
.0 1. 



on ill 



where {Jx)2 is a 2 x 2 matrix. 
We define 



(5.35) 



^2{z) = ei{z)-e2{z), 



z G 



so that we have by (5.22) 



■ gni/'2,+(2) 
(5.36) (Jx)2(z) = ( ^ gnV2,-(^) ] ' 



z G iM, \z\ > y*{a). 



Then ■ip2 is associated with the measure cj2, since 



az az 



.da2 



t(si -(z) - si,+ (2;)) = 2TTi — {z), 

az 
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by (2.7), and also 



We also identify an external field V2 on the imaginary axis, which is only 



f 



there in case a > 0. We have by (|5.23) 
(5.37) {Jx)2{z) - 

with 



g-nV2{z) ]_ 



z e {-iy*{a),iy*{a)), 



(5.38) 



Vo(z) 



62{z) - ei{z), for z G (-iy*(a), %*(«)), 
0, 



elsewhere on zM. 



The external field V2 will not be active, since it acts only on the part of ^2 
that is in (— iy*(a), iy*(Q)), and this part is zero, since S{ii2) = S{a2) = 
i'K\{-iy*{a),iy*{a)). 

See Figures [5^ and [5^ for the jump matrices Jx in the two cases a > 
and Q < 0. 
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6 Second transformation X U 
6.1 Definition of second transformation 

The second transformation of the RH problem uses the functions that come 
from the vector equihbrium problem. It is possible to state the transforma- 
tion in terms of either the 5- functions, or the A- functions. 

Definition 6.1. We define the 4x4 matrix valued function U by 



(6.1) 



Uiz) 



I + ^nCE2A ) e"^X(z)e-"^(^) 



where C is the constant from (4.31), G is given by 



(6.2) 



G = diag {gi + h 52 - 9i 93 - 92 -93) 
= - diag (Ai -V X2- 01 A3 - 6*2 A4 



3) 



and L is a constant diagonal matrix 

(6.3) L = diag {ii O) , 

with ii the variational constant in the Euler-Lagrange condition on fii . 



Note that the equality of the two diagonal matrices in (6.2) follows from 
the definition (4.35) of the A- functions. 

Then U is defined and analytic in C \ (M U iM) . 

6.2 Asymptotic behavior of U 
Lemma 6.2. We have 



(6.4) C/(z) = (/ + 0(z-i/3)) 



/I 













^1/3 














1 





Vo 








^-1/3) 



1 

A, 



as 2: —7- 00 in the jth quadrant. 

Proof. We have because of the asymptotic behavior of the A functions that 



^-nG{z) ^ ^jj^g ^^_n^-nh ^n/3 ^n/3 ^n/3^ 

3 
2 



Qj 
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as z — >■ oo in the jth quadrant, where 

' f^i = diag (l CO o;^) 

= diag (w 1 w^) 

ils = diag (w^ 1 a;) 

_ Jl4 = diag (l oj^ ui) 

Then by the asymptotic behavior of X, and the definition of U, we get 



2 lO 



as z ^ oo in the j th quadrant. We can move the 0{z-^) to the front, but 
then the ©(z'^/^) reduces to 0{z-'^/^): 

U{z)=(^I+^nCE2,^ + 0{z-y^)Jdiag{l zV3 i z'^^) (j^ ^) 



We also want to move the z term to the left. Then we pick up an 
0{1) contribution in the (2, 4) entry. Indeed we have 



z-2/3diag(l ^V3 1 ^-1/3) 





Aj^ljAj 



ijdiag(l Z-V3 1 ^1/3) 

= E2,4 + ©(^-'Z^') 



as z — >■ DO in the jth quadrant, since 



AjQjAj'^ = I 1 



r 

1 
1 0> 



for every j. The lemma follows. 



□ 



Remark 6.3. As a consistency check, we compute the jumps of the matrix 
valued function A{z) defined by 

^(z) = diag(l zV3 1 z-V3)(^l 
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for z in the jth quadrant. Then we have by the definition (5.13)-(5.14) of 
the matrices A 



J' 



(6.5) A+iz) = A-{z)JA{z), 

with jump matrix 



(6.6) 



Ja{z) 











o\ 





1 

















1 







-1 


i 


: 








o\ 








-1 








1 
















V 



z G 



which is indeed what we expect. 
6.3 Jump matrices for U 

The jump [/+ = U-Jjj with jump matrix Jjj takes a different form on the 
various parts of the real and imaginary axis. Since G and L are diagonal 
matrices, the jump matrix Jjj has the same block structure as Jx- In terms 
of the A functions the jumps take on a very nice form. 

The jump matrix Jjj on the real line has the block form 



(6.7) 



Ju{x) 



{Ju)i (x) 
(Juhix) 



X G 



with 2x2 blocks {Ju)i and {Ju)3- On the imaginary axis it takes the form 
(6.8) Ju{z) = I (Ju)o(x) 1, ze 



'1 0> 

(Jt/)2(x) 
.0 1> 



with a 2 X 2 block {Ju)2- 
Lemma 6.4. We have 



(6.9) 



{Ju)i 



gra(Ai,+-Ai,_) 





g'^(A2,+ -^l.-) 

gn(A2,+-A2,_) 



on 
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(I 

g«(->»2. + -A2.-) 

1 
\Q 







jn(A3, + -A3,_) 








1/ 



IC2 



/e"(Ai,+-Ai,_) 




o"(A2.+ -A2,-) 



-V 



/I 









0\ 




e"(^2-^3) 1 



V 



1/ 



, g-2ji7rMfc gn(A2-Ai) 



e"('*'^' + ^'^3'-^ 






1 



V 



Figure 6.1: Jump matrices Ju in case a > 



(6.10) 

and 

(6.11) 



(^1 



gn(A2,+-A2,_) Q 
gra(A2,+-A3,_) gn(A3,+-A3,_) 



,n{A3,+-A3,_) n(A4,+-A3,_) 



on i& 







on(A4,+-A4,_) 1 ' 



on 



Proof. By (|6.1|), \Q.2\ we get 



g-?i(Ai_-V') Q 

Q g-n(A2,^-6»i, 



gn(Ai,+-y) Q 

) ) (-^^^i ( e"(^2,+-ei,+) 



Now ^1 is analytic across M \ {0}, and Vi = V — Oi, so that by (5.29) we 
indeed obtain (6.9). 

For {Ju)2 find in a similar way 
(6.12) 

^g-n(A2,--9i,-) Q \ /gn(A2,+-ei,+) q 

Q g-n(A3,--e2,_) ) (•^^)2 ( Q gn(A3,+-e2.+) 



For z e iM with \z\ > y*{a) we have Oi^±{z) = 92,^{z) by (2.10) and then 
KM, and Kl^ give us that KM holds on iM \ (-y*(a), iy*(a)). 
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/I 








0\ 





gn(A2, + -A2,-) 











1 


gn(A3, + -A3._) 





Vo 








V 




^g«(Ai, + -Ai._) I Q 

e"(^2-+"^2-) 
e"(^3.+-A3,-) I 

e"(-^4,+-A4,-) 



e"-{^3,+ ->'3,-) I 



n(A4_ + -A4._) 



Figure 6.2: Jump matrices Jjj in case a < 



On {—iy* {a),iy* (a)) (which is only relevant in case a > 0) we also ob- 
tain ( |6.10| ), but now we use the fact that 9i and 02 are both analytic on 
{-iy*{a),iy*{a)), together with (IS^TJ), (fS^Sl), and ([6l2|). 



For (J[/)3 we obtain from (6.1) and (6.2) 
(6.13) 

/e-"{^3,--e2,-) 



^^3 



gn{A3,+-e2, 







pn(A4,+-e3,+) 



For (Jx)3 we have the two expressions (5.30) and (5.32). Using this in (6.13) 



we obtain (6.11) in both cases. 



□ 



The expressions in (6.9)-(6.11 ) are valid over the full (real or imaginary) 
axis. Observe in particular that the two expressions (5.36) and 
both lead to 



6.10 



and the two expressions (5.30) and (5.32 




(Jx) 

(Jx)3 both lead to (6.11). Hence the special roles that ±x*{a) (in case 
a < 0) and ±iy*{a) (in case a > 0) played in the jump matrix Jx for X 
have disappeared in the jump matrix Jjj for U. 



6.4 RH problem for U 



We have found the following RH problem for U. 
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(6.14) 



U is analytic in C \ (R U iR), 
U+ = U-Ju, onMUiR, 
?7(z) = (/ + 0(z-i/3))diagfl 1 z-l 



1 
A, 

as z — 7- oo in the j'th quadrant, 



where Ju is given by (6.7)-(6.8) with {Ju)k A: = 1,2,3 given by (6.9)- 
The parts {Ju)k ™ jump matrix Jjj have different expressions in the 



various parts of the real and imaginary axis. This follows from ( 6.9 )-( 6.11) 



and the jump properties of the A-functions as given in Lemma 4.12 Also 
recall that n is a multiple of three. 



Lemma 6.5. (a) For {Ju)i we have 
(6.15) 

''e"{-^i,+-^i,-) 1 ^ 

g"{A2,+ -A2,-) 

g -2n7riofe gn(A2,+-Ai,_)\ 



{Ju)i 







on S{fii), 

on {bk,ak+i), 
fork = 0,l,...,N. 



(b) For {Jjj) 2 we have 
(6.16) 

^gn(A2,+-A2,_) Q ^ 

I g"(A3,+ -A3,-) 

gn(A2,+-A3 _) Y j ' 



U)2 



on S{a2 - /U2), 

on iR \ S{(T2 — /i2), 



(c) For {Ju)3 we have 



3.17) {Ju)3 = { 



gn(A3,+-A3,-) I 

Q g"(A4,+-A4,-) 
g»^{A4,+-A3,_)\ 

(-1)" J ' 



, on 8(^13), 
onR\S{fi3). 



Proof. All expressions follow from (6.9)-(6.11) and Lemma 4.12 
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In particular, note that on (—03,03) (which can only be non-empty if 
a < 0), we have by (|4.42) 



„n(A3,+-A3,_) _ n(A4,+-A4,-) _ „n(-|7ri) 



i-iy 



since n is a multiple of three. This explains the entries (—1)" in (6.17) on 

(-C3,C3). □ 



7 Opening of lenses 

The next step in the steepest descent analysis is the opening of lenses around 
S{fii), S{a2 — fJ-2) and S{fi3). This will be done in the third and fourth 
transformations U ^ T and T ^ S. 



7.1 Third transformation U t-^ T 

In the third transformation we open the lenses around S{a2 — ^2) and S{fi3) 
that will be denoted by L2 and L3, respectively. The lenses will be closed 
unbounded sets that do not intersect. There are three situations, depending 
on whether C2 and C3 are positive or zero. We recall that by regularity we 
can not have 02 = 03 = 0. The three different cases differ in the shapes of 
the lens, which is due to the different supports of S{a2 — 1^2) and 8(113) as 



illustrated in Figures 7.1 7.2 and 7.3 



In the Cases IV and V we have that both C2 and C3 are positive and 
we choose the lips of the lenses such that they have ±102 and ±03 as 



endpoints, as shown in Figure 7.1 



In Case III we have C2 = and C3 > 0, and now we open the lens such 
that the lips around S{a2 — /^2) = stay away from the imaginary 



axis and intersecting the real line at two points ±73 as in Figure 7.2 
We recall that in Case III we have ^ S{fii) and hence the number 
N of intervals in S{fj,i) is even. We then choose 73 such that 

(7.1) < 73 < min(a^/2+i) C3), in Case III. 

The lens around S{p3) is as in the Cases IV and V. 

In the cases I and II we have C2 > and C3 = 0. We then take the lens 
around 5(o"2 — /i2) as in the Cases IV and V above, but we choose the 
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Figure 7.1: Unbounded lenses in case C2 > and C3 > 0, corresponding to 
Cases IV and V. The figure also shows the jump contour and some of 



the jump matrices in the RH problem (7.9) for T. 



lips of the lens around S^fi^) = M such that it is away from the real 
axis and intersects the imaginary axis at two points ±472 with 

(7.2) < 72 < C2, in Cases I and II. 

We choose all lenses to be symmetric with respect to both the real and 
imaginary axes. 

Note that we have 73 in Case III and 72 in Cases I and II. For ease of 
presentation we also define 

72 = in Cases III, IV, and V, 

(7.3) 

73 = in Cases I, II, IV, and V. 

There are further requirements on the lenses that are important for the 
steepest descent analysis. These are formulated in the next two lemmas. 
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Figure 7.2: Unbounded lenses in case C2 = and C3 > 0, corresponding to 
Case III. The figure also shows the jump contour and some of the jump 
matrices in the RH problem (7.9) for T. 



Lemma 7.1. We can and do choose the lens L2 around S{a2 — ^2) such 
that 

Re(A3 - A2) < in L2 \ iR 

and such that 

{z £ C \ \z\ > R, \ Re{z)\ < e\ Im(z)|} C L2 
for some e > and R > 0. 



Proof. From (4.35) we obtain 

A3 - A2 = 2g2 -91-93 + ^2-01 
and for z G S{a2 — ^2) we have by ( 4.11[ ) and (2.10) that 

(A3 - \2)±{z) = ± (<72,+ (^) - 52,-(^) + eiA^) - 9i,+ {z)) 
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Using (4.9), (2.7) and (2.9) we can further rewrite this to 



(A3-A2)±(z) 



zbivrz lb 27ri((T2 — /^2)([0, -z]) for Imz > 0, 



=FoVri =F 2vrz((T2 — /i2)([-Z5 0]) for Imz < 0, 



which is purely imaginary. Then 

ay \dz\ 

which is positive for y G M in Case III, and for \y\ > C2 > in the other 
cases. 

The first statement of the lemma then follows from the Cauchy Riemann 
equations since Re(A3 — A2)± = on the imaginary axis. The second state- 
ment follows then from the asymptotic behavior of A2 and A3 as given in 
(|4.47l), (|4.48l) and Lemma [214} □ 



Lemma 7.2. We can and do choose the lens L3 around S{iJL-i) such that 

Re(A3 - A4) < in L3 \ M 

and such that 

{z e€.\\z\> R,\ Im(z)| < e| Re(z)|} C L3 
for some e > and R > 0. 



Proof. The proof is similar to the proof of Lemma 7.1 



Using (4.35) we obtain 

A3 - A4 = -2^3 + 52 + 6*2 - O3. 

and hence 

A3,± - A4,± = -53,+ - 53,- + g2,± + 02,± - 93,± ± (53 



93,- 



Now by (2.15) and (5.33) 



V3{x) ±27ria3{[x*{a),x]), x > 0, 
V3{x)^2ma3{[x,-x*{a)]), x < 0. 



Combining this with (4.12) and ( |4.7[ ) leads to 

0-3 ([x* (a), x]) - ^3([x,oo)), 



A3,±(3;) — A4,±(x) = ±27ri x 



X > C3, 



-cr3([x, -x*(a)) - ;U3([x, 00)) - |, X < -C3, 
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Figure 7.3: Unbounded lenses in case C2 > and C3 = 0, corresponding to 
Cases I and II. The figure also shows the jump contour Sj- and some of the 
jump matrices in the RH problem (7.9) for T. 



which is purely imaginary. We find 

|:Im(±(A3-A4k(x)) = 2.^fe±i^(x) 

and this is positive for x G M in Cases I and II, and for |x| > C3 > in the 
other cases. 



As in Lemma 7.1 the first statement of the lemma now follows by the 
Cauchy Riemann equations. The second statement follows by (4.48), (4.49), 
and Lemma \2M □ 



Now that we have defined the lenses L2 and L3 we can come to the actual 
definition of the transformation U ^ T. The transformation is based on the 
following factorization of {Ju)2 on S{(J2 —^2) as given by (6.16) (recall that 
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(7.4) 



A3,t) 



(Juh 



g?l(A3-A2)- 

1 

]_ gn(A3-A2)- 

1 





pn(A3-A2)H 

-1 

1 



I gri(A3-A2) + 
1 



and the factorization of {Ju)3 on ^(/xs) as given by (6.17) (recall that A3, 
A4^zp or As^-t = A4_zp lb Ittz and n is a multiple of three) 



± 



U)3 



(7.5) 



j7l(A3-A4)^ 



1 



1 

n(A3-A4)_ 

1 
-1 



1 

gn(A3-A4)+ ^ 



g?l(A3-A4)^ 

This leads to the following definition of T. 

Definition 7.3. We define the 4x4 matrix valued function T by 

r/i 



(7.6) 



T = U X { 



(7.7) 



U X < 



1 





1 
1 


\0 

/l 

1 





1 
1 




o\ 

_gn(A3-A2) Q 

1 

w 

0\ 

gn(A3-A2) Q 





V 






1 

on(A3-A4) 











1 








yO e"(^-^"^4) I J 



in the left part 
of the lens 
around S{a2 — fJ-2), 

in the right part 
of the lens 
around S{a2 — ^12)- 



in the upper part 
of the lens 
around 8(^13), 

in the lower part 
of the lens 
around 8(^13), 



and 
(7.8) 



U elsewhere. 
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Then T is defined and analytic in C \ S-p where is the contour con- 
sisting of the real and imaginary axes and the lips of the lenses around 
S{a2 - ^2) and 5(/i3). 



7.2 RH problem for T 

Now T solves the following RH problem 



(7.9) 



T is analytic in C \ S-p, 
r+ = T^Jt, on T,T, 
T{z) = (/ + 0(z-i/3))diag(l 



1 

Z3 



1 



_ 1 
Z 3 



1 

A, 



as z — )■ 00 in the jth quadrant, 



with certain jump matrices Jt that will be described in the next subsection. 



The matrices Aj are given in (5.13)-(5.14). The contour and some of 



the jump matrices Jt are shown in Figures |7.2| and |7.3 

Since the lenses around S{a2 — ^2) and S^j.^) are unbounded, we have to 
be careful about the asymptotic behavior of T(z) as z — )• 00, since it could 
be different from the asymptotic behavior of U{z). However, a simple check 



using (4.47)-(4.49 ), (7.6) and (7.7) shows that the asymptotic behavior is 



the same. 

However, it is good to note the following. The asymptotic behavior of the 



Pearcey functions given in (5.9) is not uniform up to the real and imaginary 



axes. By following the transformations Y ^ X ^ U we see that the same 
is true for the asymptotic behavior oiU . It requires an independent check 
that after opening of the unbounded lenses the asymptotic behavior of T is 
in fact uniform in each of the quadrants. This phenomenon also appeared 



m 



7.3 Jump matrices for T 

Our next task is to compute the jump matrices Jt in the Riemann-Hilbert 
problem (7.9) for T. The definitions (7.6) and (7.7) of T and the structure 
of the jump matrix Ju as given in (6.7) and (6.8) yield that 



T+ = T-Jt 
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where Jt has again the structure 



(7.10) 



Jt 



{Jt)i 
. (Jt)3, 




on M \ (-73,73) and 
on the hps of the lens 
around S^fis), 

on iR. \ (— Z72, ^72) and 
on the hps of the lens 
around S{a2 — IJ-2), 



with 2x2 blocks {JT)k for k = 1,2,3. The block structure as in (7.10) is 
not valid on the intervals (—73, 73) and (—^72, ^72) (if non-empty). On these 
intervals the block structure changes to 



(7.11) 




on (-73,73), 
on (-^72,^72), 



with some non-zero entries that are denoted by *, see (7.15) and (7.16) 
below. 

The diagonal blocks are given in the next lemma. 



Lemma 7.4. (a) For (Jt)i we have 
(7.12) [Jt)i - 



{Ju)i on M, 

I2 on the lips of the lens around S{p,3), 



where {Ju)i is given by (6.15). 
(b) For (Jt)2 we have 



(7.13) (Jt)2 = < 



\ ^ri{\3-\2) 

,0 1 ^ 

1 



2,+ - 



-^3,-) 1 



on 5(0-2 - M2), 

on the lips of the lens 
around S{a2 — /U2), 

on {-ic2, ic2). 
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(c) For (Jt)3 uie have 



(7.14) (Jt)3 = < 




on S{fj.3), 

on the lips of the lens 
around S^fis), 

on (-C3,C3). 



Proof. The jumps follow from the jumps Ju for U in (6.15 )-(6.17) and the 



transformation (7.6)-(7.8), where we use the factorizations (7.4) and (7.5). 

Outside the lenses the blocks (Jt)2 and ( Jt)3 have not changed, and so 
on {-ic2,ic2) = iM. \ S{a2 - ^12) we have ( Jr)2 = {Ju )2 an d on (-C3, C3) = 
M\ S{iJ,3) we have (Jt)3 = {Ju)3- See also ( |7.16| ) and (7.15) for the result of 



the calculation of the full jump matrix Jt on the intervals (—^72,^72) and 
(-73,73)- □ 



It remains to describe the off-diagonal entries in (7.11) that only occur 
in case 73 > or 72 > 0, that is in Cases III or I/II. 

Case III: 73 > 0. In Case III we have C2 = and the lips of the lens 
around S{a2 — ^2) intersect the real line at ±73 with 73 > 0. 

In order to compute the jump matrix Jt on (—73,73), we have to note 
that by the regularity assumption we have that ^ S {fj,i)U S (fis) , and by the 
choice of 73 in (7.1) we have that [—73,73] is disjoint from ^(/xi) H 8(113). 
Also, the measure fii is symmetric so that ajv/2 = 1/2- Hence the jump 



matrix Jjj given in (6.7), (6.15) and (6.17) takes the form 



Ju 



/(- 




gn(A2,+-Ai,-) 








\ 







(-ir 


















(-ir 


gn(A4,+— As,- 


) 


V 











(-1)" 


/ 



on 



-73,73)- 



Prom (7.6) we obtain 



Jt 



/I 0\ 

1 zpe"(^3-A2) 

1 

Vo 1/ 



Ju 



/I 0\ 

1 ±e"(^3-A2) 

1 

Vo 1/ 



on 



-73,73) n : 
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After some calculations using (4.37), (4.38), (4.42) and the fact that n 



is a multiple of 3 we obtain from the last two expressions that 



(7.15) Jt 



g"{A2,+-Ai,_) _|_gn(A3,+-Ai,_) 

(-1)" 

(-1)" 

V 



\ 

n(A4,+-A2,-) 
A3,-) 

(-1)" / 

on (-73,73) n 



=Fe 

n{A4 



si 



Cases I and II: 72 > 0. In Cases I and II we have C3 = and the lips of 
the lens around S{fj,3) intersect the imaginary axis at ±272 with 72 > 0. By 



(6.8) and (6.16) we have 



J, 



u 



/I 



e 

Vo 



0\ 

1 

"(A2, + -A3,-) I Q 

1/ 



on (-i72,i72)- 



Combining this with (7.7) and using (4.40), the fact that A4 is analytic on 



\ {0} and the fact that n is a multiple of 3 leads to 



(7.16) 



Jt 



(\ 0\ 

1 

e"(^2,+-A3,-) 1 Q 

Vo ±e"(^2,+-A4,-) \j 



on 



-i72, ^72) n ill 



7.4 Fourth transformation T ^ S 

In the fourth transformation we open up a lens around each interval [ofc, h]^ 
of S{fii). The union of these lenses will be denoted by Li and Li is a closed 
bounded set. This is done in a standard way based on the factorization 



(Jr)i = (Juh 



(7.17) 



gn(Ai-A2)+ I 

Q gn(Ai— A2) 

1 

gn(Ai-A2)- ^ 



1 
-1 



1 0^ 

gn(Ai-A2)+ Y 



Lemma 7.5. We can and do choose the lenses Li around S{ni) such that 

Re(Ai - A2) < in Li \ M. 
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Figure 7.4: The opening of the lenses around S{fii). The lens around (a^, bk) 
intersects the lens around S{fi3) if and only if ±03 € {ak,bk) as is the case 
in the figure. 



Proof. By (4.35) we obtain 

(Ai - A2)± = V -h- 25i,± + c/2,± - Oi. 



By (4.10), (2.14) and this leads to 



\^{gi,+ {x) - gi,^{x)) = ±27rifii{[x,oo)), x > 0, 

I =f(5i,+ (x) — lb ovri = ib27ri/ii([x, cxd)) it ovri, x < 0. 



Thus ±(Ai — A2)± is purely imaginary and 
d 



which is strictly positive on each of the intervals {ak,bk) by the regularity 
assumption on fii. The lemma follows by the Cauchy Riemann equations. 

□ 
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In addition to the condition described in the lemma we also make sure 
that the lenses around S{fii) do not intersect the lens around S{a2 — H2) and 
the lips of the lens around ^(/ia), except for the case when ±03 G S{fii). In 
that case we choose the lips around the interval(s) containing ±03 in such 
a way that they intersect the lips of the lens around S{fi3) exactly once in 



each quadrant as it is shown in Figure 7.4 



If € S{fii) then the lips of the lens around the interval containing 
intersect the imaginary axis in the two points where 

< 72 < C2. 



See also Figure 7.4 If in addition G S{fis) then also 

< 72 < 72. 

If <S'(/ii) then we put 

72 = 0. 

We also take the lenses so that they are symmetric in both the real and 
imaginary axis. 

We now define S as follows. 

Definition 7.6. We define 

' ( 1 o\ 

_en{Ai-A2) 10 

10 





(7.18) 



S = Tx I 



1 



1 
o\ 



I V 



en(Ai-A2) 10 

10 








V 



in the upper part 
of the lenses 
around S{^i), 

in the lower part 
of the lenses 
around S{fii), 



and 



(7.19) 



S 



elsewhere. 



Then S is defined and analytic in C\T,s where S5 is the contour consist- 
ing of the real and imaginary axis and the lips of the lenses around S{fj,i), 
S{a2 - /U2) and 5(^3). 
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7.5 RH problem for S 

The asymptotic behavior as z — )• oo 
satisfies the following RH problem. 



clearly has not changed, and so S 



(7.20) 



S is analytic in C \ S5, 
S+ = S-Js, on S5, 
S{z) = (/ + 0(z-i/3))diagfl 



1 



1 



_ 1 

Z 3 



1 
A, 



as z ^ 00 in the jth quadrant 



with jump matrices Js that are described next. 

The jump matrices Js have again the block structure 



(7.21) Js = { 




/l 0\ 

{Js)2 

I Vo 1/ 



on M \ (-73,73) and 
on the lips of the lenses 
around S{fj,i) and S{fxs), 

on iR \ ((-^72, ^72) U (-■ 172, ^72)) 
and on the lips of the lens 
around S{a2 — IJ-2), 



with 2x2 blocks {Js)k for k = 1,2,3. The block structure is different on 
the intervals (—73,73), (— ^72,^72) and (— i72,i72) if non-empty. 
On these intervals we have 



(7.22) 




on (-73,73), 



on (-Z72,n2) or {-1^2,1%), 



with possible non-zero entries that are denoted by *. 
The diagonal blocks are given in the following lemma. 
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Lemma 7.7. 

(7.23) 



(a) For {Js)i we have 



s)i 



-1 



1 






1 

^n{Xi-X2)(z) 



(b) For {Js)2 we have 



(7.24) {Jsh 




(c) For {Js)3 we have 
(7.25) 



(J: 



S)3 




h 



1 







3?l(A3-A4) 



on S{ni), 

on {bk,ak+i) 

fork = 0,...,N, 

on the lips of the lenses 

around S'(^i), 

on the lips of the lens 

around S{fi3). 



on S{a2 - ^2), 

on the lips of the lens 
around S{a2 — ^2), 

on {-ic2,ic2). 



on 5(^3), 

on (-C3,C3), 

on the lips of the lenses 
around S{fii), 
on the lips of the lens 
around S{fis). 



Proof. The expressions for {Js)i on S{fii) and on the hps of the lenses 
around S{^i) fohow from the factorization ( |7.17 ) and the transformation 
(7.18)-( [7.19 ). On the other contom's we have {Js)i = {Jt)i- Also {Js)k = 
{JT)k for ^ = 2, 3, and so we find all further expressions in the lemma from 
([7!Il|-((7l4l). □ 
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We next give the jump matrices Js on the intervals (—73, 73), (— ^72, ^72) 
and (—^72,^72) and in particular the off-diagonal entries that were denoted 



by * in (7.22). They depend on the Cases I-V. We simply present the 



formulas without further comment. Of course they follow from the jump 



matrices Jt and the tranformation (7.18)-(7.19). Note that 



72 > <^ G S{iJ,i), 73 > 44^ C2 = 0, 72 > 44> C3 = 0. 



Case I. In Case I we have G S^fii), C2 > and C3 = so that 72 > 
and 72 > and 73 = 0. By construction 72 < 72 and we find in the Case I 



(7.26) 



Js={ 



( / 1 o\ 

10 

-pgn(Ai-A3_) gn(A2,+-A3,_) ^ q 

^ _gn(Ai-A4) _|_gra(A2,+-A4) q 

[Jt 



on (-272,^72) n il 



on (-^72, U (i72,'i72), 



where Jt is given by (7.16). 



Case II. In Case II we have S{fii), C2 > and C3 = so that 72 = 0, 
72 > and 73 = 0. In this case 



(7.27) 



Js = Jt on(-i72,i72) 



where Jt is given by (7.16). 



Case III. In Case III we have 5'(/xi), C2 = and C3 > so that 72 = 0, 
72 = and 73 > 0. We have in this case 



(7.28) 



Js = Jt on (-73,73) 



where Jt is given by (7.15). 



Case IV. In Case IV we have E S{fii), C2 > and C3 > so that 72 > 0, 
72 = and 73 = 0. We have in this case 



(7.29) Js 



/I 0\ 

10 

-pg7i(Ai-A3,_) gn(A2,+-A3,_) Q 

V 1/ 



on (-^72,^72) n il 
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Case V. In Case V we have S{fj.i), C2 > and C3 > so that 72 = 0, 
72 = and 73 = 0. There are no exceptional intervals in this case. 

7.6 Behavior of jumps as n — ?► 00 

Having collected all the jump matrices Js we may study their behavior as 
n — )• 00. It turns out that all off-diagonal entries of the form ±e"^'*'^~^'=^ are 
such that Re(Aj — Afc) < and therefore they are exponentially decaying as 
n — )• 00. 

For {Js)i in (7.23) we have off diagonal entries 



(7.30) 



g"(A2,+-Ai,-) 
gn(Ai-A2) 



on M\5(/xi), 

on the lips of the lenses around S{fii), 



which are indeed exponentially decaying because of (4.43) and Lemma 7.5 



For {Js)2 in (7.24) we have off diagonal entries 
(7.31) 



gn(A2,+-A3,_) 
gn(A3-A2) 



on iR \ S{(T2 — /U2), 

on the lips of the lens L2 around S{a2 — ^2)1 



which are exponentially decaying because of (4.44) and Lemma |7.1 
For {Js)3 in (7.25) we have off diagonal entries 



(7.32) 



gn(A4,+-A3,_) 
gn{A3-A4) 



onM\5(/U3), 

on the lips of the lens L3 around S^^s) 



which are exponentially decaying because of (4.45) and Lemma 



The remaining off-diagonal entries appear in the * entries in (7.22 ) on the 



7.2 



special intervals (—73,73) or {—1^2,1^2) and (— i72,i72)- They are explicitly 
given in the formulas (7.16), (7.15), (7.26), and ( 7.29| ). It turns out that all 
these entries are exponentially decaying as n — )• 00. We will not verify all 
the cases here, but let us check the jump matrix J5 on (—^72,^72) in Case 



I as given in (7.26). Here there are four off diagonal entries 

^n(Ai-A3,_) n(A2,+-A3,_) _ n(Ai-A4) 



(7.33) 



±e 



n{A2,+-A4) 



Since we are in Case I, the inequalities from (4.44) and Lemmas 7.2 and 7.5 



apply on [— ^72, ^72] since this interval is contained in iM \ S{a2 — /U2) and it 
belongs to the two lenses Li and L3. So we have 



(7.34) Re(Ai - A2) < 0, Re(A2,+ - A3,-) < 0, Re(A3 - A4) < 
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on [—ij2,i'%]- In fact, equality in the first and third inequalities of (7.34) 



holds only at 0. Then indeed all entries in (7.33) are exponentially decaying 
as n — >■ oo, uniformly on [—^72,^72]- 

In the next step of the steepest descent analysis we will ignore all expo- 
nentially small entries in the jump matrices Js- This will lead to matrices 
Jm that we use as jump matrices for the so-called global parametrix. The 



matrices Jm are given in (8.3)-(8.6) below. 



8 Global parametrix 



If we ignore all entries in the jump matrices Js that are exponentially small 
as n — )• 00, we find the following model Riemann-Hilbert problem for M : 

C \ (M U iM) 



■^4x4 



' M is analytic in C \ (M U iR), 

M+=M^Jm, onMUiM, 

/I \ 

M{z) = {I + 0{z-^)) I J U 

\o z-V3y 

as z — >• 00 in the jth quadrant, 



1 

A, 



where Jm is given as follows. On the real line, the jump matrix has the 
block form 



(8.2) 
with 



Jm{x) 



{Jm)i{x) 
{JmHx) 



.3) (Jm)i = < 



-1 0; 

2n7rjQfc 



on S{p,i), 











p2n7rjafe 



^■4) (Ja/)3 = < 



on (6fc,afc+i), 
for A; = l,...,iV-l, 
J2 on (-00, ai) U (^AT, +00), 

[(-l)"/2 On(-C3,C3). 
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On the imaginary axis the jump matrix has the block form 



(8.5) 



with 



(8.6) 



'1 







.0 







M)2 



0> 



Jm{x) = {Jm)2{x) 



X £ ill 



on S{a2 - fJ.2), 
on {-ic2,ic2). 



m 



Note that we have strenghtened the O term in the asymptotic condition 
1| from 0(2-1/3) to 0(^-1). 



The solution has fourth root singularities at all branch points ak,bj^, for 
k = 1, . . . ,N, and at ±C2, ±03 if 02,03 > 0. We give more details on the 
construction in the rest of this section. 



8.1 Riemann surface as an M-curve 

We follow the approach of [67J in using meromorphic differentials on the 
Riemann surface as the main ingredient in the construction of the global 
parametrix. Recall that 7^ is a four sheeted cover of the Riemann sphere. 
We use 

to denote the canonical projection, and we let ttj be its restriction to the 
jth sheet. Then for z G C we have that 

denotes the point on the jth sheet that projects onto z. It is well-defined 
for z G C\ (MUiM). 

The Riemann surface TZ has the structure of an M-curve. It has an 
anti-holomorphic involution 

(/> : 7^ ^ 7^ : P P 

with P on the same sheet as P. The fixed point set of consists of 3 + 1 
connected components {g is the genus of TZ) 

{P GTZ: (piP) = P} = So U Si U • • • U 

where 
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• So contains c«i and consists of the intervals (—00, ai] and [^atjOo) on 
the first and second sheets. 

• Sj for i = 1, . . . , N — 1 contains the intervals [bi, Oj+i] on the first and 
second sheets. 

A full description of the curves depends on the case we are in. Recall 
that there are five cases. 

Cases I and II: In Cases I and II we have C3 = and the genus is — 1. 
In these cases we have only the above curves Sq, . . . , Sjv-i and they 
are on the first and second sheets only. 

Case III: In Case III we still have g = N — 1. Since supp(//i) the 
number N is even (by symmetry). The curve Sjv/2 is now of a different 
character, since it visits all four sheets. Indeed we have that S7V/2 
consists of [6Ar/2, OAr/2+1] on the first and second sheets, and [— 03,03] 
on the third and fourth seets. The component is also unusual 

in that it contains four branch points, in contrast to the other curves 
that have two branch points. 

Cases IV and V: In Cases IV and V the genus is N. In these cases there 
is an additional component E^v that consists of the intervals [—03,03] 
on the third and fourth sheets. The other components are on the first 
and second sheets only, as in Cases I and II. 

We need the following result on non-special divisors on M-curves. 

9 

Lemma 8.1. If Pj £ Sj for j = l,...,g then the divisor ^2 non- 
special. 

Proof. This can be found in |581 Theorem 2.4]. □ 

9 

Recall that the divisor ^ Pj is non-special if and only if there are no 

i=i 

non-constant analytic functions on 7^ \ {Pi, . . . , Pg} with only simple poles 
at the points Pi , . . . , . 

8.2 Canonical homology basis 

The Riemann surface has a canonical homology basis {^1, . . . , Ag] Pi, ... , P^} 
that we choose such that the cycle Aj is homologous to Sj but disjoint from 
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<■ . 



'CLN bN 



IC2 



ai 



-IC2 



IC2 



-IC2 



Figure 8.1: Case II 



it. See Figures 8.1-8.3 for an illustration of the curves for the cycles with 
their orientation in the various cases. In all cases we choose the i?-cycles 
such that as sets they are invariant under the involution cj). 



The Cases I and II are very similar and we only show the Figure 8.1 for 
Case II. The -B-cycles are on the first sheet only. The cycle Bj surrounds 
the interval [ai,&j] once in the negative (counterclockwise) direction, and 
it intersects the real line in (— oo,ai) and in {bj,aj+i). The A-cycles are 
partly on the upper half of the first sheet and partly on the lower half of the 
second sheet. The cycle Aj passes through the cuts [oj,&j] and [aj+i,6j+i] 



with orientation as indicated in Figures 8.1 In Case II there is an A-cycle 
that intersects with the imaginary axis. We make sure that on the second 
sheet it does so in the interval {—ic2, 0). 

In Case III we have C2 = and C3 > 0. The genus of 7^ is g = — 1. 
Note that in this case, the number of intervals in S{fii) is even and the 
origin does not belong to The canonical homology basis is chosen as 
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in Figure |8.2[ The S-cycles are the same as in Case II. In particular they 
are only on the first sheet. The ^-cylces are also the same as in Case II, 
except for the cycle ^7v/2 which crosses the imaginary axis. This cycle visits 



all four sheets as indicated in Figure 8.2 



The two Cases IV and V are again very similar and we only show the 



Figure 8.3 for Case V. Here we have C2 > and cs > 0, and the Riemann 



surface TZ has genus g = N. The cycles Ai, . . . ,^Ar_i and Bi, . . . ,-BAr_i 
are as in the previous Cases I and II. There are two extra cycles Aj^ and 
Bn. The cycle A^^ is on the third and fourth sheets. It consists of a part in 
the upper half plane of TZ^ from a point in (ca, oo) to a point in (— oo, — C3), 
together with a part in the lower half plane of TZ^ that we choose to be 
the mirror image in the real line of the part on the third sheet. The cycle 
i?7V is on the second and third sheets. The part on 7^2 goes from a point in 
(— ioo, —ic2) to a point in (ic2, 00) and it goes around all intervals [aj, hj] that 
are in the left-half plane. The part on TZ^ intersect the real line somewhere 



92 



M. Duits, A.B.J. Kuijlaars, and M.Y. Mo 




-C3 C3 

^ '-"An 



Figure 8.3: Case V 

in (0,C3). 

By construction we have 

(8.7) ^{Bj) = -B„ ^{Aj)r^A,, j = l,...,g, 

where the symbol ~ means that (piAj) is homologous to Aj in TZ. 

8.3 Meromorphic differentials 

Let us now recall some facts about meromorphic differentials on the Riemann 
surface. Most of the results that we will be using can be found in [53j . 

A meromorphic differential with simple poles only is called a meromor- 
phic differential of the third kind. A meromorphic differential of the third 
kind is uniquely determined by its ^-periods and the location and residues 
at its poles, provided that the residues add up to zero. 
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We pick points Pj, j = 1,. . . ,g with Pj G Sj. For each such choice we 
define a meromorphic differential of the third kind as follows. 

We use ooi to denote the point at infinity on the first sheet and 002 to 
denote the other point at infinity which is common to all three other sheets. 

Definition 8.2. Let Pj € T,j for j = !,...,(/. Then cop is the m,eromorphic 
differential of third kind on IZ which is uniquely determined by the following 
conditions 

(a) The meromorphic differential has simple poles at aj,bj, j = 1, . . . , N , 
at ±ic2 (if C2 > 0), at ±03 (if C3 > 0), at the points Pj, j = 1, . . . ,g 
and at 002. The residues at the finite branch points are equal to —1/2: 

Res wp = Res = — i, j = l,...,N, 

z=aj z=bj 

(8.8) Res wp = -5, (only if C2 > 0), 

^ ' Z=±tC2 

Res ujp = — 2, (only if cz > 0), 

the residue at 002 is equal to 2, 

(8.9) Res ojp = 2 

Z = C02 

and the residue at the points Pj is equal to 1: 

(8.10) Resa;p = l, j = l,...,g. 



(b) The meromorphic differential has vanishing A periods in Cases I, II, 
IV, and V: 

(8.11) f Up = 0, k = 1, . . . , g, (in Cases I, II, IV, and V). 
JAk 

In Case III all A-periods are vanishing, except the one of A]\f/2- 

(8.12) ® LOp = — 7ri(5fe^jv/2) k = 1, . . . , g, (only in Case III). 



A simple count shows that the residues of wp add up to and there- 
fore the meromorphic differential up is indeed uniquely defined by the pole 
conditions and the A-period conditions. 
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If one or more of the Pj's coincide with a branch point, then the residue 
conditions (8.8) have to be modified appropriately. For example, if Pj = Uj 



then 

Reswp = h. 
p. ^ 

In this way, the meromorphic differential ujp depends continuously on the 
Pj's and is well defined for each choice of Pj G Tj, j = 1, . . . , g. 

The anti-holomorphic involution (p is used to map a meromorphic dif- 
ferential a; to a meromorphic differential (jy^{uj) in an obvious way. If u) is 
equal to fj{z)dz for a meromorphic function fj on sheet j, then 0*(a;) is 
equal to 

fj{z)dz 

on sheet j. A crucial property is that wp is invariant under 
Lemma 8.3. Let Pj G T,j, j = 1, . . . , g. Then 

Up = 4>^{ujp). 

Proof. Since all the poles of wp are invariant under the involution (p, the 
meromorphic differential (jy^{u:p) has the same poles and residues as wp. 
We have to show that their ^-periods are the same. We have 



(8.13) i (t)*{up) = f 



top. 

We have that (^[Ak) is homologous to Ak in IZ, but in the process of 
deforming ^(Afc) to A}, we pick up residue contributions from the poles of 

bJp. 

For the cycles Aj. that are only on two sheets (which is the typical situ- 
ation) we pick up a residue contribution from the two endpoints of a gap in 
the support of either or and from P^. As the deformation from (j){Ak) 
to Ak will result in clockwise loops around these points, the total residue is 
\ + \- l = Q. It follows that 



(j)*{ijjp) = j) wp = 0, 

since the ^-period oi ujp is zero for A-cycles that are only on two sheets. 

In Case III there is a cycle which is on all four sheets. If we deform 
0(^Af/2) iiito Ajqi2 we pick up residue contributions from the endpoints 
~C3, C3 and from P]^j2- Then we have residue ^ four times 
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and —1 once, so that the total residue contribution is 1 (As the residue 
contribution comes from clockwise loops around these points). It follows 
that 



(^"^{ujp) = (f> ujp + 27ri = — vrz + 27ri 



-7r«, 



since the period is — vri by definition (8.12). 



So the ^-periods of (jf^{ijOp) and oop agree, and the lemma follows. □ 
Proposition 8.4. The B-periods of up are purely imaginary and the map 

(8.14) ^' : Si X ••• X Eg ^ (M/Z)9 : 

(Pi, . . . , Pg) H> ^ ( / UJp,..., I UJp 

is a well-defined, continuous bijection. 




Proof. This follows as in [67|. The fact that the map is well-defined and 
continuous is proved as in [6T, Proposition 2.3]. Due to the fact that the 

9 

divisors Yl Pj with Pj S are non-special, see Lemma 

i=i 



8.1 



the argument 



in the proof of [67| Theorem 2.6] gives first the injectivity of ^. Then 
the invariance of domain argument of the same proof yields the surjectivity 
of □ 

8.4 Definition and properties of functions Uj 



Due to Proposition 8.4 there exists a choice of points Pj, j = I, . . . , g such 



.15) 



that the corresponding meromorphic differential up satisfies 

LOp = —2n-Kiak, k = 1, . . . , N — 1, 

d) ijjp = —mri, (only in Cases IV and V), 
Jbn 

where the equalities hold modulo liriZ. Note that up is varying with n. 



We consider n as fixed and work with up satisfying (8.15) throughout this 
section. Of course, up also satisfies the conditions given in Definition |8.2[ 

We are going to integrate up along paths on TZ that start from ooi (the 
point at infinity on the first sheet) and that on each sheet remain in the 
same quadrant. So the paths do not cross the contours Sj, j = 0, . . . , g and 
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also do not intersect the imaginary axis except along the cut S{a2 — /U2) 
that connects the second and third sheets. There may be a choice in the cut 
that one takes when passing from the first to the second sheet. However, 
this will lead to the same value for the integral because of the vanishing of 



the integral of ojp over the cycles Ak and (f){Ak), see (8.11). Note that the 



exceptional ^d-period (8.12) in Case III does not play a role here. 



With this convention for the paths we define functions Uj as follows. 

Definition 8.5. For j = 1, 2, 3, 4 and 2; G C \ (M U iM) we define Uj{z) as 
the Ahelian integral 

Uj{z) = / cop. 

J OOl 

Abusing the notation, we also write 



Uj{z) = / LOp 

J OOl 

where z is considered as a point on the jth sheet. The path from ooi to 
z £ Mj follows the convention described above, namely that on each sheet 
it stays in the same quadrant. 

Then the functions uj are defined and analytic on C \ (M U iM) with the 
following jump properties. In what follows we write = to denote equality 
up to integer multiples of 2Tri. 

Lemma 8.6. (a) For x £ (a^, bk) with k = 1, . . . , N , we have 

(8.16) ui,±{x) = U2,^{x). 

(b) For X S a^+i) with k = 0, . . . , N , we have 

(8.17) ui^+(x) = ui^_(x) — 2n7riafc, 

unless Pk is on the first sheet and x = Ti{Pk) 

(8.18) ii2,+(x) = U2-{x) + Imriak + iri, 

unless Pk is on the second sheet and x = TT{Pk) 

where 60 = —00, aj\f+i = +00, qq = 1 and ajy = 0. 

(c) For X G (—03,03) we have 



(8.19) us,+ {x) = U3^^{x) - {n + l)7ri, 

unless Pk is on the third sheet and x = TT{Pk) 
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(8.20) U4^_|_(x) = Ui-{x) + mri, 

unless Pk is on the fourth sheet and x = vr(Pfc) 

where k = N/2 in Case III and k = N in Cases IV and V. 

(d) For X G (—00, — C3) U (cs, 00) we have 

(8.21) u:i^±{x) = Ui^^{x). 

(e) On the imaginary axis we have 



(8.22) 


tii^+(x) 


= ""1 -(2;) 




/or 2; G 


^, 




(8.23) 


ti2,+ (x) 


= U2-{x) 




/or 2; G (- 




iC2), 


(8.24) 


U2,±{x) 


= ^3,t(^) 




for X G {- 


-ZOO 


-102] U [ic2,ioo) 


(8.25) 


U2.,+ {x) 


= U3_(x) 


+ TTl 


for X £ {- 




«C2), 


(8.26) 


ti4_+(x) 


= Ui-{x) 


+ vri 


for 2; G ilE 


L 





Proof. This is a straightforward but rather tedious verification. All identities 
or equivalences come down to the calculation of a period of a closed loop on 
TZ for the meromorphic differential up. 



For example, to prove (8.26) for x G iM''", we note that 



U4^+(x) — ■U4^_(x) = / UJp 

'c 



where C is the cycle on IZ that is shown in Figure |8.4| for the Cases IV 
and V. This cycle can be deformed to a sum of —(j){AN)., the cycles —Aj., 
j = 1, . . . , N — 1 and a closed loop around —102- This indeed leads to 



.26) since ujp has vanishing ^-periods and <f>{A) periods, and the only 
contribution comes from —ic2 which gives us ibvri. 

In the Case III (not shown in the figure), the cycle C can be deformed 
into the cycles, —Aj, j = 1, . . . N — 1, that include the exceptional cycle 
—Aj\f/2- It is now because of (8.12) that we obtain (8.26) for x G iM^. 

The other relations follows in a similar way. □ 

We state without proof the behavior of the functions uj near the branch 
points. They follow from residue conditions (8 

Lemma 8.7. (a) For j = 1,2 and k = 1, . . . 

Uj{z) -- 
Uj{z) ■■ 



, N, we have 
llog{z - Ok) + 0{l), asz^ak, 
llog{z -bk) + 0{1), asz^hk. 
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Figure 8.4: Cycle C 

(b) For J = 2, 3 and C2 > 0, we have 

Uj{z) = -\ \og{z ^ ic2) + 0{X), as z ^ ±ic2. 

(c) For j = 3,4 and C3 > 0, we have 

■"j W = -i log(^ T C3) + 0(1), as z^ ±C3. 

(d) If Pk is on the jth sheet of the Riemann surface, then 

Ujiz) = log(z - 7r(Pfc)) + asz^ 7r(Pfc), 

where k = 1, . . . ,g. 

(e) ^.s z — > 00, 

= C(l/z) and = -| logz+C(l), for j = 2,3,4. 

If Pfe coincides with one of the branch points then the behavior should 
be modified. This should be obvious and we do not give details here. 
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8.5 Definition and properties of functions Vj 

The functions Vj are the exponentials of the functions uj which we again 
consider as functions on C \ (R U iM.) . 

Definition 8.8. For j = 1,2,3,4 and z £ C\{RU iR) we define 

Vj{z) = e"^(^) 



The jump properties of vj follow from Lemma 8.6 We state them in a 
vector form. 

Corollary 8.9. The vector-valued function [vi V2 v^j is analytic in 
C \ (M U iM) with jump property 

[vi V2 ^4)^ = [vi V2 V3 V4) _ Jv on M U iM 

with a jump matrix Jy that takes the following form. 

(a) On the real axis the jump matrix has the block form 

with 2x2 blocks 

iJv)i = (J J) on S{ni), 

(8.27) ) /. 

for k = 0, . . . , N , and 

{Jv)3 = (J q) on Sifi-i), 

(8-28) ) \^ 

[Jv)3={ Q (-1)"/ 0^(-C3,C3). 

(b) On the imaginary axis the jump matrix has the block form 

AO 0\ 

Jv = [O {Jv)2 on iR, 

\0 -1/ 
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with the 2 X 2 block 



{Jv)2 = ( ? i ) on 5(0-2 -/i2), 



,1 0, 

^.29) 

{Jv)2=iQ on {-ic2,ic2). 



Note that each of the blocks {Jv)j, j = 1)2,3 has determinant — 1. 



From Lemma 8.7 and Definition |8.8| we obtain the behavior of the Vj 



functions near the branch points and other special points. 
Corollary 8.10. (a) For j = 1,2 and k = 1, . . . , N , we have 

Vj{z) = O ({z - aky^/^) , asz^ttk, 



Uj{z) = o{{z- bk) ^^^^ , as z^ bk- 



(b) For J = 2, 3 and C2 > 0, we have 



0({z^ic2) ^^'^] , as z ±ic2. 



Uj(Z) 

(c) For j = 3,4: and C3 > 0, we have 

Uj{z) = o{{z^ cs)"^/^) , as z^ ±C3. 

(d) If Pk is on the jth sheet of the Riemann surface, then 

Vj{z) = 0{Z - TT{Pk)), aSZ^TT{Pk), 

where k = 1, . . . ,g. 

(e) As z ^ 00, 

vi{z) = l + 0{l/z) and Vj{z) = 0{z~'^^^), for j = 2,3,4. 

8.6 The first row of M 

We now define the entries in the first row of M. 
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Definition 8.11. We define 
(8.30) Mn{z) = viiz), 

V2{z) for Imz > 0, 



1.31) Mi2{z) 

1.32) Mi3(z) 
^.33) Muiz) 



—V2{z) for Imz <0, 

vs{z) for z in the first and third quadrants, 

—vsiz) for z in the second and fourth quadrants, 

—V4{z) for Rez > 0, 

V4{z) for Re z < 0. 



It turns out that with this distribution of it signs the row vector (Mij) 
has exactly the correct jump properties that are required in the RH problem 
for M. 

Proposition 8.12. The row vector (Mij)j=i^...4 satisfies the conditions that 
are necessary for the first row of the solution of the RH problem for M . 
That is, the entries Mij are analytic in C \ (M U iM) with jump property 

(8.34) (Mil Mi2 Mi3 Mm) = (Mil Mis Mm)_Jm 



on MUiM, where Jm is the jump matrix in the RH problem for M, see (8.1), 
and also 

(8.35) (Mii(z) Mi2(z) Mi3(z) Mi4(z)) 

= (1 + 0(Z-1) 0(z-2/3) 0(z-2/3) 0(z-2/3)) 

as z ^ oo. 



Proof. This follows directly from Lemma 8.9, Corollary 8.10| (e), and Defi 



nitionETTl □ 
8.7 The other rows of M 

We will now construct the other rows of M out of the first row. 

Lemma 8.13. The vector space of meromorphic functions on TZ (including 
the constant functions) whose divisor is greater than or equal to 



i=i 



is of dimension 4. 
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Proof. Let us denote, for a positive divisor D' , the space of meromorphic 
functions on TZ whose divisor is greater than or equal to —D' by L[D'). 

Let D = 'Yfj=i Pj- ^ ^ + 3002) has a pole at 002 of order at 
most 3. Therefore, if w is a local coordinate near 002, we have 

F{w) = hnj-^ + f2W-'^ + fiw-^ + 0(1) as It; ^ 0, 

for certain numbers /i , /2 , /s • The kernel of the linear map 

L(Z) + 3oo2)^C3:F^(/i,/2,/3) 

is equal to the space L{D). Since D is non-special, the space L(D) contains 
only constant functions and therefore is of dimension 1. Therefore by the 
dimension theorem for linear maps 

dimL(L' + 8002) < dimL(L>) +3 = 4. 

However, by the Riemann-Roch theorem (see e.g. t53j), we also have 

L{D + 3002) > 4. 

This proves the lemma. □ 
We continue to use 

Let Fj be a function in L(Z) + 3oo2)- We use Fjk to denote the restriction 
of Fj to the A:th sheet and we consider the row vector 

(8.36) {MnFji Mi2Fj2 M13FJ3 MuFj^) 

Note that the possible poles of Fji^. at the points 7r(P/) are cancelled by the 
zeros of Mi^ at the same point. Therefore the row vector remains bounded 
at these points and it has the same behavior near the branch points as the 
first row. 

It is easy to check that the row vector has a jump on M U iM with jump 
matrix Jm- In order to construct the other rows we therefore aim to find 
three functions F2, F3 and F4 in L{D + 3002) such that the corresponding 
row vectors also satisfy the asymptotic condition for the respective rows 2, 3 
and 4 in the RH problem for M. This will be done in the next proposition. 



Two matrix model with quartic potential 



103 



Proposition 8.14. The RH problem for M has a unique solution that is 
constructed in the form 



M 



( Mil Mi2 Mi3 Mi4 \ 

M11F21 M12F22 M13F23 M14F24 

M11F31 M12F32 M13F33 M14F34 

VM11F41 M12F42 M13F43 MuFmJ 



Proof. Let us do this for the second row. The asymptotic condition we get 
from (8.1 ) for the second row gives us as z — >• c« in the first quadrant 

,l/3 + 0(z-2/3) 
.1/3 + 



.37) 
.38) 



Miiiz)F2i{z) 

Mi2{z)F22{z) 



v-3 



.39) 
.40) 



Mi3iz)F23{z) 
Mu{z)F24{z) 



%/3 

i 

71 



UJZ 



u;-zi/3 + 0(z-2/3) 



The first condition (8.37) is satisfied if and only if ^2(001) = 0. The 
other conditions determine the behavior of F2 near 002- 
There is an expansion 

F2{w) = fsw-^ + f2W-^ + fiW~' + 0{1) 

where w is the local coordinate near 002 which we choose to be equal to 

in the first quadrant of the second sheet. Its behavior in other quadrants 

and on other sheets is determined by analytic continuation. Thus 

(8.41) F22{z) = fsZ + /2z2/3 + /izl/3 + 0(1) 

as z — )• oo in the first quadrant. We also have 

(8.42) Mi2{z) = V2{z) = mzz-'^l^ + 7712^"^ + miz'^l^ + 0(z"^/^) 
as z — )• 00 in the first quadrant with a nonzero first coefficient 

7713 / 0. 



Then inserting (8.41) and (8.42) into the condition (8.38) we obtain a linear 
system of equations for the unknowns /3, /2, /i that has a unique solution. 
These three conditions together with the fact that i^2(ooi) = determines 
F2 uniquely. 



104 



M. Duits, A.B.J. Kuijlaars, and M.Y. Mo 



Now it requires an independent check that the conditions (8.39) and 



.40) are satisfied as well with the same function F2, and also the analogous 
conditions that come from the asymptotic condition in the other quadrants. 
This then completes the construction of the second row of M. 

The functions F3 and F4 are found in a similar way and they are used 
to construct the remaining rows of M. □ 

From the construction it follows that the entries of M and are 
uniformly bounded in n away from the branch points and infinity. More 
precisely, we have the following proposition. 

Proposition 8.15. We have that 

(a) the first columns of M and Af^* are bounded away from the points 
{oi, bi, . . . , Cat, b^}, with a bound that is uniform in n. 

(b) the second columns of M and M~* are bounded away from the points 
{oi, bi, . . . , ctv, b^}, ±ic2 and 00, with a bound that is uniform in n. 

(c) the third columns of M and M^* are bounded away from the points 
ibic2, ±C3 and 00, with a bound that is uniform in n. 

(d) the first columns of M and M^* are bounded away from the points ±03 
and 00, with a bound that is uniform in n. 

Proof. We will only proof part (a) as the others follow by similar arguments. 

Let us start with the first column of M. From the structure of the RH 
problem we see that the entries of M in the first column are analytic in 
C \ [ai, bn], bounded near 00 and have an analytic continuation across each 
interval {ak,bk) and {bk,ak+i)- Hence the entries are bounded if we stay 
away from {ai,bi, . . . ,a]\f,bj\f}. Moreover, M depends continuously on the 
parameters nuk modulo integers. By compactness of the parameter space, 
it then follows that we can choose the bound for the entries such that they 
hold uniformly in n. 

As for the first column of M~*, we note that M~* satisfies a RH problem 
that has the same structure as the RH problem for M. Then the statement 
follows from the same arguments. Alternatively, the statement for M~* 
follows from the identity 



(8.43) M- 



/l \ 

-1 

0-10 

Vo -1 J 



M, 
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where M is the solution of the RH problem (8.1) but with parameters — Qfc 
instead of for k = 1,. . . ,N. The identity (8.43) follows from the fact 
that both sides solve the same, uniquely solvable, RH problem. □ 

The following corollary will be used in the next section. 

Corollary 8.16. Let K C S{iJ,i)\{ai,bi, ... ,aiy,b]\f} be compact. Then for 
every x £ K we have 

(8.44) M-\y)M+{x)= (^^^^^^""'y^ *^ as y ^ x, 

uniformly in n. The * entries denote unimportant 2x2 blocks. 



Proof. By Proposition 8.15 the first two columns in M_|_(x) — M-|_(y) are 



uniformly bounded and since they are analytic it follows that they are of 
order 0{x — y) as y ^ x uniformly in n. Also the first two rows of Mjp^(y) 



are uniformly bounded by Proposition 8.15 Then the corollary follows, since 

M-\y)M+{x) = I + M-\y){M+{x) - M^{y)). 

□ 



9 Local parametrices and final transformation 
9.1 Local parametrices 

The global parametrix M will not be a good approximation to S near the 
branch points. Around the branch points we construct a local parametrix 
in a fairly standard way with the use of Airy functions. We will not give full 
details about the construct here but only give the relevant formulas with 
some comments. 



9.1.1 Statement of local RH problems 
Let 

BP = {flfc, 6fc I /c = 1, . . . , n} U ({±ic2, ±C3} \ {0}) . 

be the set of branch points. Note that ibic2 an ±03 are only branch points 
if they are non-zero. 

There is a possibility that ±03 coincides with one of the end points a^, b^ 
of the support of ^1. This is a case that could be handled just as well but 
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the formulas are slightly different and we prefer not to give full detail in this 
case. Thus we assume 



(9.1) 



±C3 {ak,bk \k = l,...,n}. 



We take a sufficiently small disk Dp around each of the branch points 
p € BP. The disks are mutually disjoint. We also make sure that the disks 
Daf. and 1);,^. are small enough such that they do not intersect with the lips 
of the global lens around S{fis)j ^-nd similarly the disks D±c^ do not intersect 
with the lips of the lenses around S{fj,i). Also the disks around ±102 are 
small enough so that they do not intersect with the lenses around S{fii) and 
5(^3). 

We use 

D= \J Dp 
peBP 

to denote the union of the disks. Then we would like to find a parametrix 
P : D \ S5 —7- C^^^ that has the same jumps as S on the parts of in D. 
There is a very minor complication here in case that or 6^ belongs to M\ 
SifJ-s) = (— C3, C3). Then the jump matrix J5 for S has in its right lower block 
(-^5)3 an off-diagonal entry that is exponentially small in a full neighborhood 



of Ofc or see the second formula in (7.25). Being exponentially small this 



entry will play no role in what follows. However, for the construction of the 
local parametrix it is more convenient to set this entry equal to 0. This is 
also what we did when defining the part ( Jm)3 in the jump matrix for M, 
see (8.4). Thus we use (Jm)3 for the jump matrix (Jp)3 on the real line in 
Daf. and Dfe^., see (9.6) below. 

A similar thing happens in case ±03 G M \ S{^i). Then an exponentially 



small entry is in {Js)i, see the second formula in (7.23) that we set equal to 
zero in the jump matrix for P. 

Thus in Dp where p is one of ak,bk, for A: = 1, . . . , A^, we take 



(9.2) 
with 

(9.3) 



Jp 



{Jp)i 




{Jp)i = {Js)i 

{Jp)3 = {Jm)3 
{Jp)3 = h 





{JP)3 



on E5 n Dp, 

on Rn Dp, 

on {T.s\M)nDp 
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If C2 > then in Dic2 and -D-ic2 we put 

/I 0^ 

(9.4) Jp = (Jp)2 

\0 1, 

with 

(9.5) ( Jp)2 = ( Js)2 on S5 n Dp. 



If C3 > then in Dc^ and D-c^, we have again the block structure (9.2), 
but now we put 

{Jp)i = {Jm)i onMn(Z)c3UL>_c3), 
(9.6) (Jp)i=/2 on (Ss\M)n(De3UZ)_c3), 

(Jp)3 = (^5)3 on n (Dc3 U Z)_C3), 

With this definition of Jp we look for a parametrix 

that satisfies the following local RH problem. 

P is analytic in D \ and continuous in D \ S5, 

(9-7) < P+ = P^Jp, on ^sr\D, 

P = {I + 0{n^^))M as n ^ 00, uniformly on dD \ T^s- 

9.1.2 Airy functions 

Around each branch point the construction of P is essentially a 2 x 2 problem, 
that can be solved in a standard way using Airy functions, see [33] for the 
2x2 case and for larger size RH problems. 

The RH problem for Airy functions is the following. It will be stated in 
terms of an auxiliary C, variable. 



/^-l/4 



is analytic, 
on Sa, 









(1/4 



1 



(/ + 0(r3/2)) 



2^3/2 
63^ 



as C 



00, 
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Figure 9.1: Contour and jump matrices J a in the RH problem for A. 



where the contour and the jump matrices J a are shown in Figure 9.1 
We have stated the RH problem in such a way that detA(Q — )• 1 as 
— )• oo, which implies that also the solution has constant determinant 1. 
This explains the factors \/27r and iti that appear in the solution (9.9) below. 
Define the three functions 

yo(C) = Ai(C), yi(C) = ^ Ai(wC), 2/2(0 = Ai{u^C), 

where Ai denotes the usual Airy function. That is, Ai is the unique solution 
of the Airy differential equation y"{C) = CuiC) with asymptotic behavior 



Ai(C) 



_2a3/2 

e 3'' 



1 + o{c 



-3/2^ 



as C — oo with — vr + e < arg^ < vr — e for any e > 0. Also yi and y2 are 
solutions of the Airy equation and the relation 



yo + yi + 2/2 = 



is satisfied. 



Then the solution of the RH problem (|9.8|) is given by 

for < argC < 27r/3, 
for 27r/3 < argC < it, 
for — vr < argC < — 27r/3, 
for - 27r/3 < argC < 0. 



(9.9) A{C) = V27T X <^ 



yo(C) 


-y2(C)\ 


-%o(C) 


J 


-yi(C) 


-y2{C)\ 


Vi(C) 


W'2iC) J 


-y2(C) 


yiiC) \ 






yo(C) 


yiiC) ' 


-w'oiO 


-%'i(C), 
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The solution A of the Airy RH problem is the main building block for 
the local parametrix P. 



9.1.3 Parametrix P in D 



In the neighborhood D^j, of a right-end point of the support of fii the 
local parametrix P takes the form 



(9.10) 



P{z) = M{z) 



P(z) 




h 



zGDt,\^S. 



To describe the 2x2 block Pi{z) we use the solution A{() of the RH 



problem (9.8), the functions Ai and A2 that come from (4.35) (and that also 
appear in the jump matrix (Jp)i = (Js)i, see (7.23)), and a function 



2/3 



(9.11) 



f{z) 



{liXiiz)-X2{z))±l7Tiak) 



2/3 



if bk > 0, 
if bk < 0, 



for iblmz > 0, z e Db^ \ M. 

It turns out that / has an analytic extension to Db^^ which maps bk to 
(to check this one uses the jump properties ( 4.37 )~( 4.38) for Ai and A2, 
among other things). Shrinking Dh^. if necessary, one has that C = f(z) is a 
conformal map from Df,^ to a convex neighborhood of C = with f{bk) = 0, 
f'{bk) > and which is real for real z. We adjust the lens around {ak,bk) 
in such a way that the lips of the lens within are mapped by / into the 
rays arg ( = 27r/3. 

Then Pi is given by 



(9.12) 



1 / 1 



A /ni/6/(2)i/4 



X A{n^/^f{z)) 



''g|n(Ai(^)-A2(2)) 








n-i/6/(^)"^/' 



g-|n{Ai(^)-A2(2)) J ' 

for z G Db^ \ S^, ±Imz > 0, 



with the principal branches of the fourth roots in /(z)^^/^. 
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9.1.4 Pjirametrix P in Da,^ 



The construction of P in a neighborhood Da^ of a left endpoint Ofe of S{iJ,i) 
is similar. Here we have a map 



(9.13) f{z) 



'(f(Ai(2)-A2(z))±|7riafc_i)'/' if > 0, 

. ~ ^2(^;)) ± |7riafc_i =F ^Tri)^^^ if < 0, 



for zblmz > 0, z G Z^a^ \ M, with /(a^) = 0, /'(a/j) < 0. If necessary, we 
shrink the disk Da^. and adjust the hps of the lens around (a^, bk) such that 
C = 7(2;) is a conformal map in Da^. onto a convex neighborhood of ^ = 
that maps Eg H Da^. into S^. 

Then the local parametrix takes the form 
















/g|n(Ai(z)- 


-A2(Z)) 


V 





2; G £>„fe \ S5, 



nV6/(2)V4 

n-V6/(^)-i/ 

g-in(Ai(z)-A2(z)) I 
for z G L'afe \ S5, iblmz > 0. 



The fourth roots /(-z)^^^^ are taken to be analytic in Da,^ \ [0^,6^] and 
positive for real z < ak, z E Da^. 

9.1.5 Parametrix P in -D±ic2 

This case is only relevant if C2 > and so we assume C2 > 0. The parametrix 
in D±ic2 takes the form 

/I 0^ 

(9.16) P{z) = M{z) P2{z) 

Vo 1, 



where P2{z) is a 2 x 2 matrix valued function in -D±ic2 \ ^2- 

In Dic2 it is constructed with the function 
(9.17) 

/3 1 
/(^) = ( 4(-^2(^) - A3(^)) T -vrij for z G Ac2 \ ±Re2; > 0, 
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which (with an appropriate understanding of the 3/2-power) has an exten- 
sion to a conformal map on Dic^ with /(ic2) = and f'{ic2) G iR.'^. Then 
f{z) > foT z = iy e Dic2 n iM, y < C2- We take P2 in Dic^ as 



(9.18) P2{z) = i-1 
X Ain^/'fiz)) 




„l/6j(^)l/4 

n-V6/(z)-V4 

|n(A2(z)-A3(2)) 







e2 







for z G Ac2 \ ^5- 



The fourth roots /{z)"^^/^ are positive for z = iy E Dic^ D iW, y < C2- 

The construction of P in -D-ic2 is very similar. By symmetry we can 
also obtain it from P in Djc2 by means of the formula 



P{z) = 
for z e D_ic2 n E5. 



/I 












n 













1 








Pi-z) 





1 














-1 











-1 













y 












V 



9.1.6 Parametrix P in D±c^ 

This case is only relevant if C3 > and so we assume C3 > 0. The parametrix 
P in -D±C3 takes the form 



(9.19) P(.) = M(.)(;^ ^^°^,) 

where a 2 x 2 block P^iz). 

It is constructed with the function 
(9.20) 

/3 1 
/(z) = f - {Xsiz) - X^iz)) ± -TTi j for ± Im z > 0, z G D-C3 \ M. 

Then / has an analytic extension to I?-C3 which maps ±03 to 0. Shrinking 
D-I-C3 if necessary, one has that ( = f{z) is a conformal map from D-C3 to 
a neighborhood of = with /(— C3) = and f'{—c^) > 0. We adjust the 
lens around (—00, — C3) in such a way that the lips of the lens within -D_C3 
are mapped into the rays dixg( = 27r/3. Then P3 is given by 
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(9.21) P,{z) = {-1) 



1 / 1 



V2 



A /ni/6/(^)i/4 







-i 1 J \ n-i/6/(2)-i/4^ 

g|n(A3(^;)-A4(2)) q \ 

Q g_l„(A3(^)-A4(2)) J ' 



for z e D^C3 \ ^5- 



For the construction in Dc^ we can use the conformal map / with the 
same definition as in (9.20), but now considered in the neighborhood of C3. 
Then we have /(cs) = and /'(cs) < 0, and after adjusting of the lenses we 
then define P3 as 



(9.22) P,{z) = {-1) 



X A{n^/^f{z)) 



1 /I 



V2 

1 

-1 



-i\ (n^l^f{zfl^ \ 

n-V6/(^)-i/4^ 

^ ^\n{\-i{z)-M{z)) Q \ 

f,-\n{\-i{z)-Xi{z)) I ' 



for z G \ Ss- 



9.2 Final transformation 

Having the global parametrix M and the local parametrix P we are ready 
for the fifth and also final transformation. 

Definition 9.1. We define the 4x4 matrix valued function R by 
(9.23) 

^ I (iisfo around each of the branch points, 

]S{z)M{z)^^ outside the disks. 

Then R is defined and analytic outside the union of and the bound- 
aries of the disks around the branch points. However, since the jumps of 
5* and M agree on 5(/ii) n S{fi3) and on S{a2 — ^2), the matrix R has an- 
alytic continuation across the parts of these supports that are outside the 
disks. We also have that the jumps of S and P agree inside the disks D±ic2 
and so R has analytic continuation in the interior of these two disks. In 
the disk Dai^ ,Di,^ and D±c^ the jumps of S and P may not be exactly the 
same. They could differ by an exponentially small entry on the real line 
inside these disks. The jumps on the lips of the lenses are the same inside 
the disks, and so R has an analytic continuation across these lenses, but R 
could have an exponentially small jump on the real line inside the disks. 
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The result is that R has an analytic continuation to C \ S^j for a certain 
contour S/j and that R satisfies the following RH problem. 



(9.24) 



where 



(9.25) 



R is analytic in C \ S/j, 

R+ = R-Jr, on S/^, 
^ R{z) = I + 0{z~^) as z 



oo. 



Jr{z) 



M{z)P{z)-^, 
P^{z)Js{z)P+{z)-\ 



M_(z)Js(z)M+(z)-i, z G SfiVUpesp^P- 
All of the jump matrices are close to the identity matrix as n 



oo. 



Indeed because of the matching condition in (9.7) we have 

Jr{z) = Ii + 0{n~^) as n — > oo 

uniformly for z G Upesp dDp. The other jumps are exponentially close to 
the identity matrix as n — )• oo with a bound that improves as z — )• oo. Indeed 
we have 

Jr{z) = I4 + C'(exp(— cn(|z| + 1)) as n ^ 00 



uniformly for z G T,r \ (UpeBP^-^pj- This follows as in [l5l section 8.4]. 



Observe that from the definition (9.23) and the asymptotic behaviors of 



S and M as given in ( |7.20[ ) and we first find that R{z) = h + ©(z^^/^) 
as z — ?■ 00. However Jp — I4 is also exponentially decaying as z — )• 00 with 



n fixed, and therefore the better bound R{z) = I4 + 0{z~^) in (9.24) indeed 
holds. 

From this we conclude as in 



Proposition 9.2. There is a constant C > such that for every n, 

\\R{z)-h\\< ^ 



n{\z\ + 1) 



uniformly for z £ C \ S/j. 



This concludes the steepest descent analysis of the RH problem for Y. 



114 



M. Duits, A.B.J. Kuijlaars, and M.Y. Mo 



9.3 Proof of Theorem lOl 



We use the following consequence of Corollary 8.16 and Proposition 9.2 
Lemma 9.3. For every x G S{fii) \ UfcLi i^a,, U -Df,j.) we have 
,-U„.\c _ lh + 0{x-y) 



S+ {y)S+{x) 



as y 



uniformly in n. The * entries denote unimportant 2x2 blocks. 



Proof. By (9.23) and the fact that x is outside the disks, we have 

(9.26) S^\y)S+{x) = M:,\y)Rl\y)R^{x)M+{x), 

if y is close enough to x. Note that R may not be analytic in a neighborhood 
of X. However, by deforming contours into the lower half plane we see that 
R does have an analytic continuation to a neighborhood of x where the same 
estimate of Proposition |9.2| is valid. Now write 

Rl\y)R+[x) = 1 + Rl\y) {R+{x) - R+{y)) . 

Then by Cauchy's Theorem 



R:,\y)R+ix) = 1 + R-+\y)^. i R{z 

X 



\z—x\=r 

y) 



2TTi 



I 1 

z — X z — y 
R(z) - I 



dz 



\z-x\=r {z-x){z-y) 



dz, 



for some r > 0. Combining this with Proposition |9.2| leads to 

'x - y" 



R^\y)R+{x) = I + O 



n 



as y 



X. 



Now inserting this in (9.26) and using Corollary 8.16 gives the statement. □ 
We follow the effect of the transformations Y ^ X ^ U ^ T ^ S ^ R 

in) i k 

on the kernel K^-^ given by (1.37). 

The transformation Y ^ X as given by (5.17) then gives 

1 



2TTi{x 



y) 



(O tL'0,n(y) W\,n[y) «^2,n(y)) 



(9.27) 







DnP-iiy)e 



1 

e-"®+(^)p^_+(x)L'-i 



n0+{y) X^'^{y)X+{x) 



/1\ 
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By (5.26) we then get 
(9.28) 



2TTi{x — y) 



(0 e-"(^(?')-^i(2')) Q)X-^{y)X+{x) 



n\ 






Then by the transformation X given in (6.1) we get 



K^^{x,y) 



(9.29) 



g"(9i,+W+^l) 

27ri(x — y) 
X U^\y)U+{x 



(O e^"(^(2/)-fi(iy)-si,+ (»/)+sr2,+ (2/)) Q o) 



/1\ 






By (4.35) this leads to 



, ^ Q'n.{g\M^)-9i,+{y)) , , 

K\l>{x,y)= — ^^—^ — (O e"(^2,+(y)-Ai,+(y)) q q) 



(9.30) 



27ri(x — y) 



X ^-ny)f/+(x) 



/1\ 







The transformation U ^ T given by (7.6)-(7.8) only acts on the lower 



right 2x2 block, and does not affect the expression (9.30) for the correlation 
kernel. We get 



2TTi{x — y) 



(O e"'^-^2,+{?/)-Ai,+(i/)) Q 



(9.31) 



xT-\y)T+{x) 



/1\ 





voy 
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Next, by the transformation T i— )• 5 of (7.18) and (7.19) we find 
(9.32) 

gn(gi,+(x)-gi,+ (j/)) 



27ri{x — y) 



X S^\y)S+{x) 



(-X5(Mi)(y) e"(^^.+(j^)-^^.+(^)) 0) 



/ 1 \ 



V / 



where Xs{^ll) denotes the characteristic fmiction of the set S{fii). 

Let X G S{fii) \ {ak,bk \ k = 1,...,N}. The factor e"(9i.+(^)-9i.+(2')) 
disappears as y ^ x. Then we use Lemma |9.3| and we get 



")(x,x) = lim + 0(1) 

(9.33) y^"" 27r^(x - y) 



n d 
2TTi dx 



(Ai,+ (x)-A2,+ (x)) + 0(l) 



as n — )• oo. 
Thus 



lim -K'f^{x,x) = 77^;^ - A2,+ (x)) 

n^oo n Im dx 



1 



(Fi,_(x)-Fi,+ (x)) 



where we used (4.35), (4.10), (4.13), and (4.14). This completes the proof 
of Theorem 11.41 
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